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ABSTRACT 


A  computerized  system  for  reducing  middle  atmos¬ 
phere  electrical  conductivity  data  obtained  by  blunt  probes 
and  Gerdien  condensers  is  developed.  In  addition  to  im¬ 
proving  the  computational  speed  and  accuracy,  this  particu¬ 
lar  system  enhances  the  instrument's  sensitivity  range, 
i.e.  the  altitude  region  over  which  measurements  can  be 
obtained . 


The  reduction  system  is  designed  strongly  toward 
software  portability  and  micro-processor  adaptability. 
Electrical  conductivity  data  from  recently  launched  blunt 
probes  were  reduced  using  this  system,  and  the  results  are 
discussed  in  this  thesis.  Future  expansion,  improvements 
and  two  proposed  micro-processor  system  designs  are  also 
presented . 
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CHAPTER  I 


INTRODUCTION 


J_.J_  BACKGROUND 

The  reduction  of  blunt  probe  and  Gerdien  con¬ 
denser  data  from  strip  charts  can  be  a  relatively  tedious 
and  time  consuming  task.  Also,  maunal  reduction  tech¬ 
niques  are  limited  in  accuracy.  With  the  advancements  in 
computer  technology  and  real-time  programming  techniques, 
this  particular  task  is  readily  adaptable  for  a 
computer-assisted  reduction  method.  In  fact,  such  a 
method  was  initially  attempted  at  The  University  of  Texas 
at  El  Paso  a  few  years  ago  with  modest  success  [Shih 
(1977)].  Recently,  this  research  effort  was  continued 
using  state-of-the-art  techniques  which  resulted  in  an 
operator- inter  acting  ,  computerized  reduction  procedure 
for  obtaining  electrical  conductivity  from  blunt  probe 
and  Gerdien  condenser  measurments.  The  presentation  of 
this  research  is  the  subject  of  this  thesis. 

U2  STATEMENT  OF  THE.  PROBLEM  AND  THESIS  OUTLINE 

The  computerized  data  reduction  system  incor¬ 
porating  hardware  design,  software  interfacing,  programs 
and  design  strategy  is  presented  in  this  thesis.  The  new 
system  also  provides  data  base  management  techniques  for 
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In  Chapter  2,  a  brief  description  of  the  theory 
of  operation  for  the  blunt  probe  and  the  Gerdien  con¬ 
denser  and  the  equations  for  calculating  electrical  con¬ 
ductivity  are  presented.  A  survey  of  former  data  reduc¬ 
tion  methods  is  considered  in  Chapter  3.  Chapter  4  will 
discuss  the  hardware  design  for  the  data  reduction  system 
while  Chapters  5,  6  and  7  describe  the  philosophy,  design 
strategy,  functionality  and  implementation  of  the  system 
software.  In  Chapter  8,  proposed  micro-processor  reduc¬ 
tion  systems  are  considered,  one  at  the  system  level  and 
the  other  at  the  component  level.  Finally,  Chapter  9  is 
concerned  with  the  error  analysis  of  the  data  reduction 
system,  and  Chapter  10  presents  results  and  considers  fu¬ 
ture  improvements  to  this  new  system. 


CHAPTER  II 


BLUNT  PROBE  AND  GERDIEN  CONDENSER  EXPERIMENTS 

The  blunt  probe  [Hale  and  Hoult  (1965);  Hale 
(1967);  Hale,  Hoult  and  Baker  (1968)]  and  Gerdien  con¬ 
denser  [Pedersen  (1964);  Rose  and  Widdel  (1972);  Conley 
(1974);  Croskey,  Hale  and  Leiden  (1977);  Mitchell,  Sagar 
and  Olsen  (1977)]  are  instruments  which  measure  electri¬ 
cal  conductivity  in  the  middle  atmosphere.  In  addition, 
the  Gerdien  condenser  measures  ion  mobility  and  charge 
number  density.  The  payload  is  launched  to  a  nominal  al¬ 
titude  of  75  km  using  a  meteorological  rocket.  Electri¬ 
cal  parameters  are  measured  after  the  payload  separates 
from  the  rocket  at  apogee  and  is  decending  on  a  stabil¬ 
ized  parachute.  The  probe's  current- voltage  response 
characteristics  are  telemetered  to  a  ground  receiving 
station  and  recorded  on  magnetic  tape  for  later  data 
reduction.  A  brief  description  of  the  instruments' 
operation  is  presented  in  the  following  sections. 

2.1  CURRENT-VOLTAGE  RESPONSE  CHARACTERISTICS 

The  charged  particle  current  collected  by  the 
probe's  electrode  is  described  by  Ohm's  law.  For  the 
blunt  probe,  the  collector  geometry  is  a  disk.  The  equa¬ 
tion  for  the  charged  particle  collection  current  is: 


9  ± 
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(2.1) 

The  collection  voltage  waveform  is  a  ramp  which  initially 
biases  the  collector  at  a  negative  potential  with  respect 
to  the  plasma,  and  then  sweeps  the  collector  to  a  posi¬ 
tive  potential.  A  complete  sweep  voltage  waveform  has  a 
period  of  approxmately  8  seconds  and  nominal  voltage  lim¬ 
its  of  5  V.  In  actual  practice,  the  atmospheric  electr¬ 
ical  conductivity  values  are  obtained  by  evaluating  the 
derivatives  of  the  probe's  current- voltage  response,  thus 
avoiding  the  problem  of  having  to  estimate  probe  poten¬ 
tial: 
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The  Gerdien  condenser  collector  geometry  consists 
of  two  concentric  cylindrical  electrodes.  The  collection 
voltage  waveform,  which  biases  the  inner  collector  with 
respect  to  the  outer  return  electrode,  is  a  ramp  voltage 
similar  to  the  one  described  for  the  blunt  probe.  The 
resulting  current  response  in  the  linear  region  of  opera¬ 
tion  is  described  by  the  equation: 


(2.3) 


The  linear  region  of  operation  for  a  Gerdien  con¬ 
denser  occurs  at  sufficiently  low  probe  voltages  such 
that  no  ion  mobility  group  is  completely  swept  out  of  the 
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air  sample  as  it  flows  through  the  instrument.  At  larger 
probe  voltages  where  one  or  more  ion  mobility  groups  are 
completely  collected  from  the  air  sample,  the  current 
response  displays  additional  structure  that  is  used  to 
calculate  ion  mobility  and  number  density.  The  reduction 
procedure  for  determining  these  parameters  is  discussed 
by  Pedersen  (1964),  Croskey  (1976),  Sagar  (1977)  and 
Domagalski  (1979). 

As  was  the  case  for  the  blunt  probe,  reduction  of 
the  Gerdien  condenser  data  to  determine  electrical  con¬ 
ductivity  actually  involves  measuring  of  the  derivative 
of  the  probe  current  with  respect  to  voltage: 

Z'l 

2.2  DATA  FORMAT 

Prior  to  launch,  the  probe's  current  response 
with  a  known  resistance  value  inserted  between  the  col¬ 
lector  and  return  electrode  is  received  through  the 
telemetry  system  and  recorded  on  magnetic  tape.  This 
pre-flight  calibration  current  response  to  the  known  col¬ 
lection  voltage  waveform  is  later  scaled  and  used  in 
reducing  the  in-flight  conductivity  data.  The  in-flight 
data  are  also  telemetered  back  to  ground,  preferably  to 
the  same  receiving  station,  and  recorded  on  magnetic 
tape . 


d\± 
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The  probe's  modulation  scheme  involves  converting 
the  electrometer's  analog  output  signal  to  a  pulse  train 
which  has  a  frequency  proportional  to  the  measured 
current.  The  relatively  low  frequency  pulse  train  (nomi¬ 
nally  0  -  200  Hz)  modulates  the  grid  of  the  1680  MHz 
transmitter.  This  particular  modulation  scheme  is  compa¬ 
tible  with  the  Meteorological  Rocket  Network's  GMD-3  re¬ 
ceiving  system  and  TMQ-5  strip  chart  display  system.  In 
playing  back  the  recorded  data  for  reduction,  the  first 
step  involves  continuous  measurment  of  the  telemetered 
pulse  frequency. 


For  this  particular  modulation  scheme,  the  equa¬ 
tions  for  conductivity  as  functions  of  pulse  frequency 
are  written  for  the  blunt  probe  and  Gerdien  condenser  as 
follows : 


R  I  (</f/</t)DATA 
- — - { BLUNT  PROBE ) 

2  r*  Rcal  (dVd\)  CAL 

(2.5) 


a  ± 


In  (r0/r 8 )  _ I _  {d1/dt)oATA 

2.W  X  Rcal  (di/di) cal 


(  GEOIEN  CONOENER) 


(2.6) 

In  these  expressions,  the  pulse  frequencies  ^CAL  an<* 
f D^r a  correspond  to  the  telemetered  waveforms  received 
during  preflight  calibration  and  while  the  instrument  is 
in  flight,  respectively. 
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CHAPTER  III 

DATA  REDUCTION  METHODS 

From  the  equations  developed  in  Chapter  2,  it  is 
seen  that  reducing  the  blunt  probe  and  Gerdien  condenser 
data  requires  measuring  both  (df/dt)DATA  and  (df/dt)CAL. 
Various  methods  have  been  used  to  obtain  these  values. 
Two  former  methods  which  will  be  discussed  are  the  strip 
chart  method  and  the  computerized  technique.  In  addi¬ 
tion,  the  new  computerized  data  reduction  method  will  be 
introduced  . 

3-1  STRIP  CHART  METHOD 

The  strip  chart  method  is  probably  the  most  used 
method  for  conductivity  data  reduction.  The  recorded 
waveforms  are  played  back  through  a  frequency-to-voltage 
transducer  from  which  an  analog  signal  is  obtained  and 
displayed  on  a  strip  chart.  The  waveform's  slopes  are 
then  scaled  and  measured.  From  these  measurments,  the 
conductivities  can  readily  be  computed  using,  in  addi¬ 
tion,  the  physical  parameters  of  the  probe  and  the  cali¬ 
bration  waveform  measurments.  Examples  of  data  waveforms 
displayed  on  strip  charts  are  shown  in  Figure  3.1. 

In  spite  of  the  simplicity  of  this  method,  there 
are  inherent  disadvantages  associated  with  it,  primarily, 
because  the  strip  chart  recorder  is  a  mechanical  device. 
The  pen's  time  response  degrades  the  display  and  the 
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scaling  of  waveforms  obtained  early  in  the  flight  is  very 
difficult  since  the  conductivity  values,  i.e.  the 
waveform's  slopes,  are  relatively  large.  At  the  other 
extreme,  data  obtained  late  in  the  flight  have  slopes  so 
flat  that  scaling  the  waveforms  introduces  error.  Final¬ 
ly,  noise  during  the  flight  is  usually  exaggerated  by  the 
pen's  inertia.  The  important  advantage  to  this  method  is 
that  a  hard  copy  of  the  reduced  waveforms  does  result. 

3.2  FORMER  COMPUTERIZED  DATA  REDUCTION  METHOD 

An  earlier  attempt  was  made  to  reduce  conductivi¬ 
ty  data  using  a  PDP  11/10  mini-computer  with  a  Laboratory 
Peripheral  System  (LPS)[  Laboratory  Peripheral  System 
User's  Manual  (1973)3.  Computerized  techniques  enhance 
computational  speed  and  accuracy,  and  they  increase  the 
altitude  region  over  which  the  conductivity  values  can  be 
extracted  from  the  data  waveforms.  In  the  initial  ap¬ 
proach,  however,  there  also  were  limitations. 

First,  the  package  was  not  designed  with  a  uni¬ 
form  data  structure.  Thus,  the  data  for  each  flight  were 
considered  individually  and  formatted  differently. 
Secondly,  the  timing  reference  which  was  generated  inter¬ 
nally1  would  at  times  became  distorted  if  signal  dropout 
occurred  for  extended  time  intervals.  Lastly,  large  por- 

•internal  timing  means  the  programmable  real-time  clock 
is  used  for  both  frequency  measurment  and  time  accmula- 
tion . 
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tions  of  the  programs  were  written  in  assembler  language, 
thus  making  it  difficult  to  maintain  or  to  modify.  This 
also  reduced  protability  between  different  machines. 

3.3  IMPROVED  COMPUTERIZED  DATA  REDUCTION  METHOD 

This  research  in  computerized  data  reduction 
techniques  not  only  improves  the  former  system  but  in  ad¬ 
dition,  it  extends  the  capabilities  with  respect  to  por¬ 
tability  of  software,  hardware  design  and  data  base 
management.  The  resultant  product  is  a  complete  package 
that  includes  hardware  design  and  software  program 
modules  to  perform  data  acquisition,  data  reduction,  data 
display  and  data  base  management. 

3.3-1  11SI1  PRINCIPLE 

The  basic  concept  of  computerized  data  reduction 
is  very  simple.  If  the  probe’s  transmitted  signal  ,  i.e. 
the  waveforms,  can  be  digitized,  then  they  can  be  manipu¬ 
lated  and  processed  by  computers  in  a  form  suitable  for 
reduction.  Furthermore,  they  can  be  stored,  retrieved 
and  managed  through  mass  storage  media. 

The  method  for  digitizing  waveforms  is  now  ex¬ 
plained.  The  incoming  signal,  which  is  in  the  form  of 
variable  frequency  pulses,  is  compared  against  a  rela¬ 
tively  higher  fixed  frequency  signal  (clock).  Each  time 
a  pulse  is  detected,  the  count  from  the  fixed  frequency 


source  is  accumulated  and  referenced  to  the  previously 
detected  pulse.  This  count,  which  is  a  direct  measure¬ 
ment  of  the  time  elapsed  between  two  successive  pulses, 
is  then  used  to  determine  pulse  frequency  (Figure  3.2). 
This  method  is  superior  to  using  a  combination  of  a 
frequency-to- vol tage  converter  and  an  A/D  converter  to 
achieve  the  digitization. 

After  the  waveforms  have  been  digitized,  they  are 
stored  on  disk  for  further  processing.  Each  waveform  is 
then  properly  identified  and  accompanied  by  a  header  that 
contains  relevant  information  including  physical  parame¬ 
ters  of  the  instrument,  the  time  and  date  of  the  launch 
and  the  timing  of  that  waveform.  These  waveforms  are 
written  out  to  the  disk  as  a  file  and  stored  for  future 
processing . 

The  next  step  is  scaling  the  waveforms.  Instead 
of  displaying  the  waveforms  on  a  strip  chart,  the 
waveforms  are  shown  on  a  storage  scope.  The  operator 
identifies  different  slopes  on  the  waveform  and  directs 
the  computer  to  use  numerical  methods  to  calculate  the 
best-fit  straight  line.  The  straight  line  that  is  calcu¬ 
lated  is  generated  on  the  screen  to  demonstrate  that  the 
fit  is  satisfactory.  If  the  operator  feels  that  the 
straight  line  is  the  proper  one,  the  slope  of  the  line  is 
used  for  the  calculation  and  the  results  are  stored  back 
in  the  header  of  each  waveform.  From  this  point  on,  the 
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flight  consists  of  individual,  self-contained  waveforms, 
each  complete  with  all  information  necessary. 

Further  processing  includes  plotting  and  printing 
of  the  data,  data  library  creation  and  management,  and 
transportation  of  the  data  to  other  computer  systems. 

3.3.2  MAJOR  CHANGES 

Three  major  changes  have  been  included  with  the 
new  computerized  system  for  reducing  electrical  conduc¬ 
tivity  data.  They  are:  inclusion  of  an  external  time 
code  generator  for  referencing;  implementation  of  a  fixed 
data  structure  for  computational  flexibility;  and  pro¬ 
gramming  in  higher-level  language  for  software  portabili¬ 
ty.  A  discussion  of  these  improvements  is  considered  in 
the  following  sections. 

3.3.2. 1  EXTERNAL  TIME  CODE  GENERATOR 

The  data  reduction  system  uses  an  external  time 
code  generator  for  referencing  the  data  waveforms.  The 
programmable  real-time  clock  of  the  LPS  is  thus  dedicated 
to  measuring  pulse  frequency. 

When  the  payload  launch  is  detected  on  the  data 
tape  playback,  the  time  code  generator  is  maunally  ini¬ 
tiated  by  the  operator.  For  each  data  waveform,  the 
operator  issues  a  command  via  one  of  the  A/D  converter2 
2A/D  converter  is  analog- to-digital  converter. 
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channels  notifying  the  system  to  acquire  timing  informa¬ 
tion  from  the  time  code  generator.  This  is  done  by  send¬ 
ing  hand-shaking  signals  via  the  Digital-input  and 
Digital-output  ports  of  the  LPS  through  tri-state  octal 
latches^.  The  particular  time  code  generator  in  use 
sends  out  timing  information  in  BCD1*  digits  which  are 
converted  into  binary  form. 

3. 3.2.2  FIXED  DATA  STRUCTURE 

A  fixed  data  structure  facilitates  communication 
between  program  modules  as  well  as  transportation  between 
different  machines.  Furthermore,  it  facilitates  inter¬ 
facing  and  future  modification  of  the  existing  program 
modul es . 


Each  waveform  is  accompanied  by  a  128  16-bit  word 
header  that  contains  all  the  physical  parameters  of  the 
instrument  as  well  as  information  about  the  flight.  Con¬ 
ductivity  values  are  included  in  the  header  after  each 
waveform  has  been  reduced. 

3.3. 2. 3  HIGHER-LEVEL  LANGUAGE  PROGRAMMING 

Most  of  the  program  modules  are  written  in  the 
higher-level  language  FORTRAN.  This  facilitates  debug¬ 
ging  and  modification  of  the  software.  FORTRAN  was  also 

3Tri-state  devices  can  exist  in  an  on,  off  and  high- 
impedence  state. 

^BCD  is  Binary  Coded  Decimal. 
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chosen  because  of  its  popularity  and  I/O  flexibility. 
The  programs  are  coded  to  conform  to  STANDARD  FORTRAN 
[ PDP  11  FORTRAN  Language  Reference  Manual  (1974)].  Ex¬ 
tensions  to  FORTRAN  in  the  compiler  are  kept  to  a 
minimum . 

Yet,  some  modules  in  the  software  package  have  to 
be  written  in  assembler  language,  as  in  the  cases  of  the 
interrupt  service  routine  and  real  time  processing,  to 
quicken  computation  speed.  However,  the  use  of  assembler 
language  is  kept  to  a  minimum. 


CHAPTER  IV 


HARDWARE  SET-UP 

The  basic  hardware  for  the  data  reduction  system 
is  discussed.  The  system's  equipment  is  classified  into 
four  categories,  each  of  which  is  described  in  this 
chapter.  A  complete  hardware  set-up  schematic  is  shown 
in  Figure  4.1. 

4. 1_  COMPUTER  SYSTEM 

A  processor  with  memory  is  required  to  execute 
programs,  perform  calculations  and  serve  as  temporary 
storage.  Since  the  data  reduction  system  works  in  an  in¬ 
teractive  mode,  a  console  terminal  is  also  necessary.  A 
mini-computer  or  even  a  micro- processor  will  satisfy  this 
requirement.  In  our  particular  system,  a  Digital  Equip¬ 
ment  Corporation  (DEC)  PDPTM  11/10  mini-computer  with  16K 
words^  of  core  memory,  a  RK05  disk  cartridge  and  UNIBUS^ 
interfacing  is  used. 

t 

4.1.1  PDP  H/10  MINI-COMPUTER 

The  PDP  11/10  [PDP  11/04/05/10/35/40/45  Processor 
Handbook  (1975)]  is  a  relatively  small  computer  in  the 
DEC  PDP  11  series  of  computers,  but  it  shares  the  same 
architecture  as  the  rest  of  its  family  (Figure  4.2).  The 
processor  is  connected  to  the  peripheral  equipment  ( in- 

51  Kword  =  1,024  words. 
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Hardware  Configuration  for  the  Data  Reduction  System 
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eluding  memories)  through  a  common  set  of  wires  called 
the  UNIBUS.  At  any  time,  either  one  of  the  peripherals 
or  the  processor  will  have  control  over  the  UNIBUS  (mas¬ 
ter)  and  will  have  access  over  any  other  members  on  the 
bus  (slave).  This  kind  of  architecture  makes  interfacing 
of  the  processor  to  other  equipment  very  convenient. 

Another  feature  of  the  PDP  11  family  of  computers 
is  that  they  provide  priority  vectored  interrupt  capabil¬ 
ity.  Thus,  external  devices  can  interrupt  the  processor 
according  to  their  pre-assigned  priority,  and  control  is 
then  passed  to  a  fixed  location  in  memory  called  the 
Vector  Address  which  is  assigned  to  each  type  of  device. 
As  will  be  presented  in  later  chapters,  the  interrupt 
mechanism  is  the  basis  of  the  data  acquisition  technique 
;  the  ease  in  UNIBUS  interfacing  thus  makes  the  hardware 
simple . 

U  .J_.2  CORE  MEMORY 

I 

i  '  Memory  is  required  in  any  computer  system  to  al- 

low  the  storage  of  instructions  and  data.  In  our  system, 
.  I  16K  words  (each  16  bits)  of  core  memory  provide  non- 

I  ti 

volatile  storage  for  the  system. 

4.K3  MASS  STORAGE  DEVICE 

Input  data  eventually  must  be  stored  in  some  kind 
of  mass  storage  device.  Disks,  drums,  diskettes,  or 

i  , 
i 

tf 

I 


tapes  can  serve  this  purpose.  One  consideration  is  the 
data  transfer  rate  of  the  mass  storage  device  which  has 
to  be  compatible  with  the  data  acquisition  rate.  If  not, 
some  incoming  data  will  be  lost.  Our  system  is  equipped 
with  a  RK05  disk  cartridge  that  can  provide  up  to  a  180K 
bytes/second  data  transfer  rate  and  2.5  million  bytes  of 
storage  capacity  [PDP  11  Peripherals  Handbook  (1975)]. 

M.2  WAVEFORM  PROCESSING,  MEASUREMENT  AND  DISPLAY 
HARDWARE 


The  clock,  counter,  Schmitt-tr igger ,  A/D  convert¬ 
er  and  display  unit  (Tektronix  603  storage  scope)  are  the 
basic  components  for  performing  frequency  counting  and 
data  display.  The  clock  generates  a  fixed  frequency 
reference  for  counting  the  frequencies  of  the  incoming 
pulses.  The  counter  accumulates  the  counts  from  the 
clock  until  a  Schmitt-tr igger  is  fired  by  the  incoming 
data  pulse  and  the  count  is  then  transfered  to  memory. 
The  counter  is,  of  course,  cleared  after  the  data 
transfer  and  it  is  then  made  ready  for  the  next  pulse. 

The  A/D  converters  serve  two  purposes.  First,  if 
the  input  data  are  fed  to  a  frequency-to-voltage  convert¬ 
er  which  in  turn  is  connected  to  an  A/D  converter,  a 
crude  representation  of  the  waveform  can  be  displayed  on 
a  storage  scope  for  waveform  monitoring.  Secondly,  the 


A/D  converters  provide  input  ports  for  operator  commands, 
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which  can  direct  the  program(s)  to  perform  selective 
functions . 

During  the  data  acquisition  phase,  the  display 
unit  indicates  the  presence  of  an  incoming  signal.  In 
the  data  reduction  phase,  the  display  unit  displays  the 
waveform  images  and  allows  the  operator  to  identify  and 
select  proper  slopes  of  the  waveform  for  electrical  con¬ 
ductivity  calculations. 

The  clock,  Schmitt-tr igger ,  counter  and  A/D  con¬ 
verters  can  be  built  using  standard  SSI  and  MSI^  com¬ 
ponents.  A  storage  scope  or  raster  scan  can  serve  as  the 
display  unit.  In  our  system,  however,  a  DEC  Laboraory 
Peripheral  System  (LPS)  and  a  storage  scope  are  used  to 
allow  maximum  design  flexibility. 

4.2.J_  LABORATORY  PERIPHERAL  SYSTEM  (LPS) 

The  LPS  [Laboratory  Peripheral  System  User's 
Manual  (1973)]  is  a  piece  of  0EM7  equipment  designed  by 
Digital  Equipment  Corporation.  It  consists  of  a  pro¬ 
grammable  real-time  clock,  eight  A/D  converter  channels, 
two  Schmitt- tr igger s ,  a  display  driver  that  can  be 
directly  connected  to  storage  scopes,  two  digital  I/O 
ports  (one  for  the  input,  one  for  the  output)  and  UNIBUS 
interface  control  logic  for  direct  interfacing  to  the 

^MSI  and  SSI  are  mediumTnd  small  scale  integrating  cir¬ 
cuit  respectively. 

?OEM  is  Original  Equipment  Manufacturer. 
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peripherals  (Figure  4.3). 

4.2.2  STORAGE  SCOPE 

The  storage  scope  (Tektronix  603)  [603  Monitor 

Instruction  Manual  (1972)]  is  connected  directly  to  the 
display  driver  of  the  LPS  and  provides  a  4096  x  4096  dot 
array  plotting  capability.  Plotting  and  erasing  are 
directly  under  software  control. 

4.3  PERIPHERAL  EQUIPMENT 

Additional  peripheral  equipment  are  also  needed 
for  this  particular  data  reduction  problem;  namely,  a 
time  code  generator,  circuitry  to  interface  the  time  code 
generator  to  the  LPS  and  a  pulse  generator. 

4.3-1  TIME  CODE  GENERATOR 

A  timing  device  is  needed  to  reference  the  time 

elapsed  from  the  launch  of  the  rocket  to  the  time  of  a 

particular  in-flight  data  waveform.  This  in  turn  enables 
the  operator  to  determine  the  instrument’s  position  and 
velocity  from  the  radar  data.  The  timing  device,  after 
started  by  the  operator  at  launch,  is  run  asynchronously, 

independent  of  the  rest  of  the  system,  to  keep  track  of 

the  time  of  flight.  An  Eldorado  1710  BCD  time  code  gen¬ 
erator  is  used  for  this  purpose  [Technical  Manual,  Model 
1706  and  1710  Time  Code  Generator  (1970)]. 
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Figure  4.3  Laboratory  Peripheral  System  Block  Diagram 
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When  a  new  waveform  is  recognized  by  the  opera¬ 
tor,  a  signal  (a  5V  rising  pulse)  is  sent  through  one  of 
the  A/D  converter  channels  to  tell  the  data  acquisition 
program  that  timing  information  is  needed.  At  the  same 
time,  an  endthenal8  is  padded  to  the  data  to  allow  a 
later  program  pass  to  segment  the  waveforms  into  indivi¬ 
dual  units.  Our  timer  provides  32  bits  of  BCD  parallel 
output  including  seconds,  minutes,  hours,  day,  month  and 
year.  Only  hours,  minutes  and  seconds  are  needed  for  our 
application.  Since  it  is  not  possible  for  the  time  code 
generator  to  input  32  bits  in  parallel  form  to  the  exist¬ 
ing  system,  an  interface  has  been  designed  to  separate 
the  signal  into  4  bytes,  8  bits  each.  (This  is  explained 
in  the  following  section.)  By  choosing  one  byte  at  a 
time,  it  is  suitable  for  micro-processor  interfacing 
since  most  second  generation  micro-processors  are  8-bit 
machines . 


The  output  of  our  time  code  generator  is  in  BCD 
form  and  therefore  must  be  decoded  into  binary  form.  It 
can  be  done  very  simply  by  software.  If  a  binary  output 
timer  is  used  instead,  decoding  can  be  eliminated. 


8An  endthenal  is  a  special  character(s)  to  indicate  the 
end  of  a  record . 
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4.3.2  TIME  CODE  GENERATOR  INTERFACE 

In  order  to  separate  the  32-bit  parallel  output 
of  the  time  code  generator  into  4  bytes  of  8  bits  each, 
four  Intel  8212  octal  latches  are  used  to  latch  the  tim¬ 
ing  information  into  the  Digital  input  port  of  the  LPS 
(Figure  4.4).  The  interface  operates  as  follows: 

When  the  operator  sees  the  beginning  of  a  waveform 
from  the  display  unit,  a  pulse  is  sent  into  a 
predescribed  A/D  converter  channel  (channel  0). 
This  pulse  tells  the  data  acquisition  program  that  a 
new  waveform  has  come  in  and  timing  information  is 
required.  Under  program  control,  a  series  of  bits 
in  the  Digital  output  port  of  the  LPS  are  set  con¬ 
secutively  (Figure  4.5);  they  serve  as  hand-shaking 
signals  for  the  time  code  generator  interface  to 
send  in  successive  8  bits  of  timing  information  in 
seconds,  minutes,  hours  and  day/year. 

Four  different  bit  positions  of  the  Digital  output  port 
are  tied  to  the  chip  enable  ( CE )  pin  on  the  four  latches. 
Therefore  each  bit  sent  to  the  Digital  output  port  en¬ 
ables  a  different  latch  to  sample  different  timing  infor¬ 
mation.  Since  the  Intel  8212  octal  latch  is  a  tri-state 
device,  all  the  latches  that  are  not  enabled  will  be  in 
the  high  impedence  state  and  look  like  open  circuits  to 
the  rest  of  the  system. 
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Figure  4.5  Hand-shaking  Signals  between  the  Time  Code 
Generator  Interface  and  the  Laboratory 
Peripheral  System 
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4.3.3  PULSE  GENERATOR 

The  primary  function  of  the  pulse  generator  is  to 
provide  a  fixed  length  pulse  for  the  A/D  converter  as  a 
command.  Since  the  data  acquisition  program  is  essen¬ 
tially  operating  in  a  polling  mode,  the  signal's  pulse 
width  must  be  controlled  to  avoid  it  from  being  carried 
around  the  loop(s)  more  than  once  and  resulting  in  an  er¬ 
ror.  A  simple  one-shot  circuit  is  used  and  the  pulse 
width  is  controlled  by  the  RC  time  constant. 

4.4  DATA  PROCESSING  EQUIPMENT 

Processing  of  the  reduced  data,  including  analyz¬ 
ing,  listing  and  displaying,  is  presently  done  using 
peripheral  equipment  on  another  computer  system.  Thus, 
such  equipment  is  not  essential  for  data  reduction. 
Equipment  used  in  this  capacity  are  a  PDP  11/45  mini¬ 
computer  [PDP  11/04/05/10/35/40/45  Processor  Handbook 
(1975)]  operating  under  the  UNIX  Time-sharing  System 
[UNIX  Programmer's  Manual  (1975)]  and  supporting  peri¬ 
pherals,  which  include  a  plotter,  a  storage  scope,  a 
line-printer,  a  9-track  industrial  standard  magnetic  tape 
unit  and  a  DECtape  unit.  The  use  of  this  second  mini¬ 
computer  system  for  data  processing  demonstrates  the  need 
for  software  portability. 


.1 


CHAPTER  V 


DATA  BASE  DEFINITION 

In  this  chapter  and  the  following  two  chapters, 
the  software  design  of  the  data  reduction  package  is  dis¬ 
cussed.  The  description  includes  data  base  definitions, 
software  design  philosophies,  implementation  methods  and 
individual  program  modules.  From  a  present  day's 
software  engineering  viewpoint,  programs  are  combinations 
of  data  structures  and  algorithms  [Wirth  (1976)]. 
Without  consistent  data  structures,  program  modules  can¬ 
not  communicate  with  each  other  efficently.  A  pre¬ 
described  data  base  definition  is  essential  for  easy  pro¬ 
gram  modification  and  tran spor tation .  In  this  chapter, 
the  data  base  definition  of  the  different  passes  during 
the  data  reducton  process  will  be  examined,  and  their 
meanings  and  implications  will  be  explained. 

5 . 1_  DATA  STRUCTURE  IN  PASS  1 

Pass  1  is  the  data  acquisition  phase.  In  this 
pass  the  incoming  signal  is  digitized  as  previously 
described.  The  digitized  data  are  stored  under  the  file 
named  "PASS  1 . TMP" .  It  consists  of  records  of  variable 
length  each  representing  a  waveform  in  digitized  form 
(Figure  5.1).  When  a  falling  pulse  is  detected  at  chan¬ 
nel  0  of  the  LPS's  A/D  converter,  a  "0"  (null  character 
in  ASCII)  is  put  into  the  current  buffer  in  memory.  This 


29 


30 


* 


UNUSE 


OATA 


ADC  READING 


OATA 


AOC  READING 


4  BYTE 
TIMING 


VARIABLE 

LENGTH 


Figure  5 . 1 

PASS  I  OATA  STRUCTURE 


128  WORD 
HEADER 


DATA 

AOC 

READING 

• 

• 

• 

• 

DATA 

ADC 

READING 

4> 

• 

• 

• 

* 

MULTIPLE  OF 
256  WOROS 


Figure  5.2 

PASS  2  DATA  STRUCTURE 


null  character  serves  as  an  endthenal  to  identify  the  be¬ 
ginning  of  a  new  waveform.  It  is  then  followed  by  4 
bytes  of  timing  information  obtained  from  the  external 
time  code  generator.  Following  these  5  bytes  are  data 
points  of  the  waveform  in  pairs.  The  first  word  is  the 
count  from  the  programmable  real-time  clock,  and  the  next 
word  is  the  reading  from  A/D  converter  channel  1.  This 
word  is  not  used  at  present  but  is  reserved  for  future 
expansion  of  the  system  when  multi-channel  transmission 
or  a  special  signal  indicating  the  probe's  zero-potential 
crossing  may  be  stored. 

Records  are  continuously  written  to  the  mass 
storage  device  until  a  falling  pulse  is  detected  at  chan¬ 
nel  2  of  the  A/D  converter.  This  signal  ends  the  data 
acquisition  phase  and  the  last  buffer  is  filled  with 
"1071118"  (ASCII  representation  of  'FI'  for  FINISH)  to 
the  next  256-word  boundary.  This  last  buffer  is  written 
out  and  the  file  is  closed. 

The  data  obtained  in  Pass  1  are  primitive 
representations  of  the  waveforms.  The  software  of  Pass  2 
properly  segments  and  refines  the  waveforms.  The  file 
"PASS  1 . TMP"  can  then  be  deleted. 
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5.2  DATA  STRUCTURE  IN  PASS  2 

The  major  function  of  Pass  2  is  to  refine  and 
re-organize  the  raw  data  obtained  from  Pass  1.  The  data 
stream  consists  consecutively  of  a  null,  followed  by  4 
bytes  of  timing  information  and  digitized  data  point 
pairs,  which  are  now  separated  on  a  waveform-by-waveform 
basis.  Each  waveform  is  accompanied  by  a  128-word  header 
and  is  rounded  to  the  nearest  256-word  block  boundary 
(Figure  5.2).  (A  256-word  block  is  the  basic  unit  used 
in  RK05  disk  cartridges.  For  different  devices,  this 
basic  unit  could  be  different.)  This  data  structure  is 
retained  throughout  the  whole  reduction  process  and 
serves  as  the  data  base  in  the  package.  The  header, 
w;,ich  contains  all  the  vital  launch  parameters,  will  be 
discussed  as  follows. 

5.2.1.  HEADER 

The  header  is  128  words  long  and  is  organized  as 
shown  in  Figure  5.3.  The  information  may  be  classified 
into  two  types:  user  supplied  and  program  supplied. 

During  Pass  2  of  the  data  processing,  the  opera¬ 
tor  is  required  to  supply  certain  information  requested 
by  the  program.  This  information  includes:  physical 
parameters  of  the  probe;  flight  information  such  as 
launch  site,  launch  date  and  launch  time;  and  probe  in¬ 
formation  such  as  the  RcAL  ^F  values  and  the  sweep 
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voltage  parameters.  This  information  is  put  into  the 
header  and  identifies  the  flight.  The  inclusion  of  this 
information  makes  each  waveform  a  self-contained  unit. 
If  part  of  the  flight  data  were  ever  lost,  the  remaining 
data  would  still  be  identifiable. 


The  second  type  of  information  in  the  header  is 
supplied  by  the  program,  or  more  precisely,  calculated  by 
the  program.  Specific  locations  in  the  header  are 
reserved  to  hold  reduction  results  from  Pass  3  of  the 
data  reduction  package.  The  results  include  electrical 
conductivities  for  both  positive  and  negative  charged 
particles.  The  header  also  reserves  room  to  hold  more 
than  one  set  of  conductivity  values,  a  possibility  for 
Gerdien  condenser  data.  Time  information  is  obtained  by 
decoding  the  BCD  time  code.  This  in  turn  enables  one  to 
determine  the  probe's  position  and  velocity  from  the  ra¬ 
dar  data. 


5.2.2  OUTPUT  FILE  FORMAT  IN  PASS  2 

After  the  data  points  have  been  segmented  into 
individual  waveforms  and  headers  are  inserted,  the  data 
points  are  rounded  to  the  nearest  256-word  block  boundary 
by  padding  the  remaining  spaces  with  "0"s  (zeros).  Two 
other  important  parameters  are  then  entered  into  the 
header;  namely,  the  number  of  data  point  pairs  and  the 
length  of  the  waveform  in  blocks.  Then  the  waveform  is 
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written  out  to  a  file  Galled  'PASS2.TMP'  and  is  ready  for 
the  next  pass. 

5.3  DATA  STRUCTURE  IN  PASS  3 

The  software  in  the  package  is  based  upon  the 
data  structure  developed  in  Pass  2.  Pass  3  does  not  at¬ 
tempt  to  change  any  of  the  data  points  but  merely 
displays  the  waveforms,  scales  the  designated  time  seg¬ 
ments,  calculates  the  conductivity  values  and  then  in¬ 
serts  them  into  the  proper  locations  in  the  header.  Oth¬ 
er  utility  programs  then  manipulate  all  the  information 
as  specified. 

In  addition  to  inserting  values  into  the  header 
of  each  waveform,  Pass  3  also  prints  the  information  into 
an  ASCII  file  and  on  the  system  console. 
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CHAPTER  VI 


SOFTWARE  DESIGN  PHILOSOPHIES 

Three  major  concerns  are  carried  throughout  the 
design  and  implementation  of  the  data  reduction  package; 
namely,  modularity,  portability  and  small  system  adapta¬ 
bility.  These  considerations  not  only  influence  the 
final  form  of  the  software  but  also  the  implementation 
methods  as  well. 

6.J_  MODULARITY 

Modularity  not  only  means  "structure  programming" 
[Dahl,  Dijkstra  and  Hoare  (1972)]  as  in  terms  of  software 
engineering,  but  it  also  means  a  separation  of  tasks, 
each  performed  by  different  phases  of  the  package. 
Although  the  programs  are  written  partly  in  FORTRAN  and 
partly  in  assembler  language,  discipline  is  enforced  in 
coding  to  produce  clean  and  readable  code.  Furthermore, 
subroutines  and  functions  are  used  freely  to  allow  pro¬ 
gram  modules  to  be  expressed  in  a  structual  manner  in 
spite  of  the  use  of  completely  unstructured  languages 
such  as  FORTRAN  and  assembler  language. 

In  addition  to  clean  coding,  the  data  reduction 
package  is  composed  of  individual  "PASSES",  each  of  which 
performs  a  specific  task  that  contributes  to  part  of  the 
reduction  process.  Designing  the  package  in  this  manner 
makes  each  phase  a  small,  manageable  task.  It  allows  fas- 
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ter  coding  and  more  importantly,  easier  debugging. 

These  passes  communicate  with  each  other  through 
the  use  of  files  (data  sets).  The  files  are  sets  of  data 
points  which  reside  in  mass  storage  devices,  such  as  disk 
cartridges.  There  are  advantages  and  also  disadvantages 
associated  with  data  file  communication.  The  disadvan¬ 
tage  is  that  it  is  much  slower  to  communicate  between 
files  in  a  mass  storage  device  than  it  is  to  put  every¬ 
thing  in  memory.  On  the  other  hand,  there  are  several 
advantages.  First  of  all,  work  space  is  greatly  expanded 
by  providing  virtual  memory^  capability  to  increase  the 
work  space.  Secondly,  tasks  need  not  be  carried  out  in  a 
continuous  manner,  thus  permitting  a  more  flexible 
schedule.  For  example,  if  all  the  data  reduction  phases 
are  not  finished  in  one  day,  the  remaining  phases  may  be 
continued  at  some  later  time  without  the  data  being  lost 
in  the  memory.  Finally,  data  in  files  can  readily  be 
transported  to  other  installations  or  computer  systems 
for  further  processing  or  displaying.  Data  files  in  mass 
storage  media  such  as  disk  cartridges  or  9-track  indus¬ 
trial  standard  magnetic  tape  make  an  attractive  and  cost 
effective  method  for  transportation. 


^Virtual  memory  means  that  disk  space  can  be  used  for 
program  or  data  space  as  if  they  are  part  of  the  memory 
available. 
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6.2  PORTABILITY 

With  the  tremendous  increase  in  software  develop¬ 
ment  cost,  software  portability  has  become  a  major  con¬ 
cern  in  software  engineering  [Kernighan  and  Planger 
(1976)].  The  hardware  is  becoming  more  sophisticated  and 
inexpensive,  and  thus  the  use  of  higher-level  language  to 
implement  software  becomes  feasible.  The  wide  existence 
of  FORTRAN  compilers  on  computer  systems  makes  FORTRAN, 
though  far  from  the  best,  the  almost  universal  program¬ 
ming  language  in  engineering.  Most  of  the  data  reduction 
package  is  written  in  ANSI  FORTRAN,  or  more  commonly 
known  as  STANDARD  FORTRAN.  The  programs  make  minimal 
usage  of  the  extension  provided  by  our  resident  compiler. 
In  cases  where  special  system  routines  are  used,  a  user- 
written  subroutine  could  readily  be  substituted. 

The  use  of  higher-level  language  in  programming 
allows  more  imagination  in  program  design,  more  under¬ 
standable  coding  and  easier  debugging  and  modification. 
In  fact,  during  the  design  stage  of  the  package,  the 
thought  that  "We  won't  use  assembler  language  unless 
necessary."  was  kept  strongly  in  mind  at  all  times. 

Furthermore,  the  programs  written  are  highly 
parameterized  [Kernighan  and  Plauger  (1976)],  i.e.  sym¬ 
bolic  names  are  used  instead  of  actual  values.  For  exam¬ 
ple,  if  "BUFFERSIZE"  is  equated  to  4096,  all  the  programs 
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will  use  the  symbol  "BUFFERSIZE"  instead  of  the  number 
4096.  Later,  if  the  buffer  size  needs  to  be  changed  to 
some  other  value,  only  the  definition  of  "BUFFERSIZE" 
needs  to  be  changed.  This  approach  permits  much  easier 
modification  and  errors  are  less  likely  to  occur. 

In  our  installation,  both  programs  and  data  have 
been  transferred  between  the  PDP  11/45,  operating  under 
the  UNIX  Time-sharing  System,  and  the  PDP  11/10  operating 
under  RT  —  11.  The  implication  is  two-fold.  First,  the 
PDP  11/45  supports  more  peripherals  and  utility  programs, 
like  plotting  routines,  and  thus  can  provide  more  attrac¬ 
tive  outputs.  Secondly,  it  demonstrates  portability 
between  different  computer  systems  and  different  operat¬ 
ing  systems.  The  results  turns  out  to  be  satisfactory. 
In  fact,  a  rough  estimation  shows  that  about  seventy-five 
percent  of  the  codes  are  portable.  The  other  twenty-five 
percent  belong  mostly  to  the  data  acquisition  phase  which 
is  machine  dependent. 

6.3  SMALL  SYSTEM  CONSIDERATION 

The  data  reduction  package  is  designed  with  small 
system  adaptability  in  mind.  This  consideration  is  re¬ 
flected  in  both  the  hardware  design  and  the  software 
design.  For  example,  two  future  micro-processor  systems 
are  presented  later  in  Chapter  8. 
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From  the  hardware  aspect,  all  the  required  com¬ 
ponents  are  standard,  easy-to-find ,  off-the-shelf  SSI  or 
MSI  integrated  circuits.  These  components  can  very  easi¬ 
ly  be  assembled  into  a  unit  that  serves  the  same  function 
as  the  LPS.  The  LPS,  of  course,  offers  much  more  flexi¬ 
bility  in  the  design  during  the  development  and  testing 
stages  but  in  the  production  stage,  a  hard-wired  unit 
could  function  just  as  well. 

The  data  acqisition  scheme  is  interrupt  driven. 
The  interrupt  feature  is  supported  both  by  mini-computers 
and  micro-processors.  Therefore,  our  data  reduction  sys¬ 
tem  can  be  incorporated  into  any  system  by  just  rewrit- 
ting  Pass  1  of  the  data  reduction  package. 

From  the  software  aspect,  two  major  concerns  have 
been  put  into  the  design  of  the  programs:  optimzed  pro¬ 
gram  space  and  data  space.  In  another  words,  the  sizes 
of  the  programs  are  minimized  rather  than  the  speed. 
Small  systems,  either  mini-computers  or  micro-processors, 
usually  have  only  a  limited  amount  of  memory  space  as 
well  as  addressability.  Accordingly,  the  programs  also 
have  to  be  small  such  that  they  do  not  take  up  too  much 
room.  This  is  first  done  by  splitting  tasks  into 
separate  passes  and  secondly,  by  carefully  optimizing 
coding,  even  at  the  expense  of  using  extra  loops  and 
subroutine  calls.  The  data  reduction  process,  however, 
is  relatively  slow  compared  to  computer  speed,  and  thus 
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not  much  is  actually  lost. 

The  other  concern  is  minimizing  data  space.  This 
is  done  by  restricting  array  size  and  buffer  size  to  lev¬ 
els  compatible  to  smaller  systems  at  the  expense  of  extra 
I/O  activities.  In  our  particular  system,  asynchronous 
I/O  is  supported  by  the  operating  system  and  therefore 
problems  of  this  nature  do  not  arise. 
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CHAPTER  VII 


INDIVIDUAL  PROGRAM  MODULE  DESCRIPTION 

In  this  chapter,  all  the  program  modules  in  the 
data  reduction  package  are  examined.  Their  respective 
functions  and  implementation  methods  are  discussed.  In 
such  cases  that  special  techniques  or  tricks  are  used, 
they  are  pointed  out  for  reference.  The  program  modules 
can  be  classified  into  two  categories;  namely,  the  dif¬ 
ferent  passes  and  the  utility  programs.  The  different 
passes  perform  data  acquisition,  data  segmentation  and 
reduction  while  the  utility  programs  do  things  such  as 
print  and  plot  results,  extract  data  and  put  the  data 
into  different  formats. 

7*1  PASS  1 

Pass  1  is  the  data  acquisition  pass.  The  incom¬ 
ing  data  are  detected  and  digitized  by  the  method 
described  in  Chapter  3.  Its  implementation  combines  pol¬ 
ling  and  interrupt  techniques.  The  polling  scheme  func¬ 
tions  essentially  as  a  big  loop  with  the  different  A/D 
converter  channels  sequentially  sampled  to  identify 
operator  commands.  The  digitization  of  data  is,  however, 
done  by  firing  of  the  Schmitt-trigger  to  interrupt  the 
processor  whenever  a  data  pulse  is  detected.  The  pro¬ 
grammable  real-time  clock  LPSKW  has  the  capability  of 
saving  the  content  of  the  counter  in  a  register  and  of 
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allowing  the  interrupt  service  routine  to  store  this 
count  in  memory.  Therefore,  data  acquisition  is  inter¬ 
rupt  driven. 

The  digitized  data  are  written  out  to  a  mass 
storage  device.  Since  I/O  activity  takes  time,  it  is  im¬ 
portant  to  keep  it  at  a  minimum  to  reduce  the  loss  of  in¬ 
coming  data.  Two  methods  can  be  used  to  accomplish  this. 
The  first  method  is  to  increase  the  buffer  size  for  in¬ 
coming  data  so  that  more  data  can  be  stored  before  each 
I/O  activity  takes  place.  However,  this  method  was  re¬ 
jected  because  the  package  is  aimed  at  micro-processor 
adaptation  and  for  such  systems,  memory  space  is  often 
limited.  The  second  approach  involves  using  a  multi¬ 
buffer  method.  Whenever  one  buffer  is  full,  the  incoming 
data  are  directed  to  another  buffer  immediately  while  the 
first  one  is  being  written  out. 

Under  the  RT -11  operating  system,  asynchronous 
I/O  is  supported  and  can  therefore  further  enhance  the 
performance  of  a  multi-buffer  scheme.  In  our  implementa¬ 
tion,  two  buffers  of  4K  words  each  are  designated  for 
this  purpose.  A  counter  and  a  pointer  are  maintained  by 
the  program  to  keep  track  of  the  number  of  points  in  the 
current  buffer  and  the  location  of  where  the  next  point 
is  to  be  placed.  When  the  current  buffer  is  full,  the 
counter  is  reset  and  the  pointer  designates  the  next 
buffer  immediately.  Simultaneously,  the  first  buffer  is 
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written  out  asynchronously  while  the  normal  data  acquisi¬ 
tion  operation  is  resumed. 

The  LPS  provides  a  digital  input  port  and  a  digi¬ 
tal  output  port  which  allow  timing  information  from  the 
time  code  generator  to  be  latched  in  under  program  con¬ 
trol  instead  of  putting  the  signals  on  the  UNIBUS  direct¬ 
ly.  If  this  feature  were  not  available,  care  would  be 
needed  in  designing  the  time  code  generator  interface  to 
avoid  damage  to  other  peripherals  on  the  UNIBUS. 

It  is  necessary  that  Pass  1  be  written  in  assem¬ 
bler  language  and  that  it  be  specifically  constructed  for 
a  designated  computer  system  since  the  interrupt  mechan¬ 
ism  and  hardware  configuration  on  various  machines 
d iffer  . 

7 . 2  PASS  2 

The  purpose  of  Pass  2  is  to  separate  the  data  ob¬ 
tained  from  Pass  1  into  individual  waveforms  and  to  in¬ 
sert  a  header  with  each  of  them.  The  I/O  for  waveform  | 

segmentation  makes  use  of  a  4K-word  input  buffer  and  a 
4K-word  output  buffer.  The  buffer  size  can  be  changed 
easily  to  fit  different  hardware  configurations. 

Because  the  output  of  the  time  code  generator  is 
in  BCD,  a  routine  to  convert  from  BCD  to  binary  numbers 
is  required  for  decoding  the  timing  information.  This 
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particular  routine  can  decode  up  to  1  hour  59  minutes  and 
59  seconds  of  timing  code.  The  binary  timing  is  also 
displayed  on  the  LED  output  of  the  LFS  to  show  the  pres¬ 
ence  of  waveforms  and  the  relative  timing.  This  pass  is 
also  written  in  assembler  language.  It  could  have  been 
written  in  FORTRAN  using  binary  read/write. 

7.3  PASS  3 

Pass  3  of  the  data  reduction  package  is  the  most 
compilicated  part  of  the  whole  system.  It  combines  a 
display  technique,  a  file  management  scheme  and  a  numeri¬ 
cal  curve  fitting  method  to  allow  the  user:  to  examine 
the  waveform;  to  expand  a  selected  portion  of  it;  to 
identify  break  points;  to  fit  the  best  slope  between 
them;  to  calculate  results;  and  to  store  the  results  back 
into  the  header. 

In  order  to  assure  adequate  memory  space,  careful 
coding  is  used  throughout  the  design  of  thi3  pass.  Due 
to  the  large  number  of  routines  in  Pass  3,  they  are 
briefly  classified  according  to  their  functions  and  dis¬ 
cussed  in  the  following  subsections. 

7.3-1  INPUT /OUTPUT  ROUTINES 

Input/Output  activities  are  needed  whenever  a 
waveform  is  brought  into  display  on  the  screen  and  then 
written  out  after  reduction  has  been  made  with  appropiate 


information  inserted  in  the  header.  In  addition,  an 
ASCII  file  is  also  created  with  the  altitude  and  conduc¬ 
tivities  for  quick  reference. 

Since  waveforms  are  consecutively  stored  in  mul¬ 
tiples  of  256-word  variable  length  records,  the  sequen¬ 
tial  I/O  approach  is  the  most  appropiate.  The  first 
block  of  each  waveform  is  read  in  and  the  waveform  length 
is  found  from  the  header.  Then  the  remaining  portion  is 
read  in  and  the  complete  waveform  is  either  displayed  or 
skipped,  as  determined  by  the  operator.  Although  sequen¬ 
tial  I/O  is  slower  than  random  access  I/O,  it  is  consid¬ 
erably  simpler.  Moreover,  the  sequential  approach  is  ap¬ 
propiate  when  magnetic  tape  is  used  as  the  mass  storage 
device  . 

The  RT  —  11  operating  system  provides  a  system 
subroutine  called  'ASSIGN'  [PDP  11  FORTRAN/RT-1 1  Extan- 
sion  Manual  (1975)3  which  allows  the  user  to  assign  logi¬ 
cal  unit  numbers  to  devices  and  files.  The  programs 
heavily  use  this  subroutine  to  communicate  to  physical 
devices  such  as  disks  and  magnetic  tapes.  In  other 
operating  systems  where  this  routine  is  not  available, 
Job  Control  Language  or  user  written  subroutines  can  be 


substituted  for  it. 
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7.3 .2  EXPANSION  AND  DISPLAY  ROUTINES 

After  each  waveform  has  been  read  into  memory,  it 
is  displayed  on  the  storage  scope  for  viewing.  The  ex¬ 
pansion  routine  and  the  display  routine  provide  an  accu¬ 
rate  representation  of  the  waveform  on  the  scope  and  at 
the  same  time,  they  provide  a  convenient  way  for  the  user 
to  identify  end  points  for  slope  calculation. 

The  expansion  and  display  routines  work  hand  in 
hand  to  allow  the  operator  to  choose  and  display  all  or 
part  of  the  waveform  on  the  screen.  The  method  and  ter¬ 
minology  used  is  confined  to  the  CORE10  standard  proposed 
by  the  SIGGRAPH  of  ACM  [Computer  Graphics  (1977)]. 

The  Tektronix  603  storage  scope  provides  4096  x 
4096  addressable  points  in  both  the  X  and  the  Y  direc¬ 
tions.  Each  point  can  be  individually  intensified  to 
make  up  lines  and  curves.  When  each  waveform  is  first 
read  in,  it  is  completely  displayed  from  0  to  200  Hz, 
which  is  the  probe's  nominal  operating  range.  Then  the 
operator  decides  if  waveform  expansion  is  needed  for  a 
more  detailed  viewing  of  selected  part(s)  (Figure  7.1). 
If  expansion  is  needed,  boundaries  in  both  the  X  and  the 
Y  directions  are  entered  by  the  operator  through  the  con¬ 
sole.  An  example  is:  "xmin  =  100,  xmax  =  300,  ymin  =  60, 

10C0RE  is  the  proposed  standard  in  Computer  Graphics 
from  the  Special  Interset  Group  in  Computer  Graphics  of 
the  Association  of  Computation  Machinary. 
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=  120"  which  would  be  to  expand  the  horizontal  seg¬ 
ment  of  the  waveform  from  the  100th  point  to  the  300th 
point,  inclusive,  and  from  60  Hz  to  120  Hz  in  the  verti¬ 
cal  direction.  In  CORE  system's  terms:  the  "world"  con¬ 
tains  the  whole  waveform;  the  "window"  is  set  to 
[xmin,xmax]  and  [ymin,ymax];  the  "screen"  is  the  4096  X 
4096  screen  surface  with  coordinates  from  0  to  4095;  the 
"viewport"  is  chosen  such  that  it  occupies  the  complete 
screen  surface  for  viewing. 

Other  efforts  have  also  been  made  to  ease  the 
reduction  process.  For  example,  both  the  horizontal  and 
vertical  axes  are  marked  in  grids  and  numbered.  The  hor¬ 
izontal  axis  is  marked  at  every  50  points  starting  from 
the  xminth  point.  The  vertical  axis  is  marked  at  every 
20  Hz  starting  from  the  yminth  Hz.  Marking  is  done  by 
drawing  dashed  lines  across  the  screen  at  the  appropiate 
locations.  These  spacings  can  be  changed  easily  to  suit 
individual  taste.  At  the  bottom  of  each  vertical  line 
and  at  the  left  of  each  horizontal  line  is  a  numeric 
string  that  indicates  the  value  of  the  line.  They  are 
either  frequencies  in  Hertz  (vertical)  or  number  of 
points  ( horizontal )  .  Character  string  generation  is  done 
by  a  standard  FORTRAN  "encode"  statement  that  converts 
from  binary  number  representations  into  ASCII  characters. 
The  ASCII  characters  are  then  generated  on  the  screen  by 
5x7  matrix  points  as  shown  in  Figure  7.2.  The  size  of 
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the  characters  can  also  easily  be  changed. 

7.3.3  STRAIGHT  LINE  FIT  AND  LINE  GENERATION 

After  each  pair  of  end  points  on  the  waveform  has 
been  chosen  by  the  operator,  the  best  fit  straight  line 
to  the  data  points  between  the  designated  end  points  is 
determined.  The  conventional  least  square  method  is  used 
to  determine  the  best  straight  line  fit.  In  order  to  ob¬ 
tain  more  accuracy  with  the  fit,  an  iteration  of  the 
least  square  method  is  made.  First,  all  the  data  points 
between  the  two  end  points  are  assigned  equal  weights  of 
1.0  and  a  first  fit  is  constructed  with  an  equation  in 
the  form 


for  which:  f  is  frequency  in  Hz;  t  is  time  in  seconds;  P 
is  the  slope  in  Hz/sec  and  Q  is  the  intercept  in  Hz.  All 
data  points  are  then  weighted  with  respect  to  their  abso¬ 
lute  distances  from  the  first  fitted  curve  (Figure  7.3). 
A  new  weighting  factor  between  0.0  and  1.0  is  assigned  to 
each  data  point  inversly  proportional  to  this  distance. 
Thus  points  farther  away  from  the  line  have  less  influ¬ 
ence  in  calculating  the  next  straight  line  fit.  A  second 
iteration  of  this  procedure  again  assigns  new  weighting 
factors  to  further  refine  the  straight  line  fit.  This 
iterative  process  results  in  a  calculation  of  the 
waveform's  slope  over  the  designated  time  segment. 
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After  the  final  values  of  P  and  Q  are  deter¬ 
mined,  the  resulting  straight  line  fitted  to  the  data 
points  is  generated  on  the  screen  for  the  operator's  ex¬ 
amination.  The  line  generation  uses  a  Digital  Differen¬ 
tial  Algorithm  [Newman  and  Sproull  (1979)]  (DDA)  for  gen¬ 
erating  a  series  of  points  to  making  up  the  line.  The 
line  generation  algorithm  also  determines  the  equation: 

f  =  Pt  +  Q  or  t  =  (f  -  Q)/P 

(7.2) 

according  to  the  ratio  in  the  X  and  the  Y  directions  of 
the  window  to  obtain  the  most  solid  line  possible. 

If  the  operator  determines  that  the  straight  line 
fit  is  not  satisfactory,  the  process  is  repeated  with 
other  end  points  on  the  waveform  designated.  In  cases 
where  the  waveform  is  so  noisy  that  a  good  fit  to  the 
data  is  difficult  to  obtain,  a  very  small  time  segment 
where  the  waveform  is  relatively  clean  can  be  selected  to 
obtain  a  fit.  The  result  is  then  examined  and  in  most 
cases,  the  fit  will  be  satisfactory. 

7.3.4  ALTITUDE  ROUTINE 

Accompanying  each  data  tape  for  each  flight  are 
usually  radar  data  that  provide  time  and  altitude  infor¬ 
mation.  A  routine  has  been  written  to  interpolate  the 
altitude  of  each  waveform  from  the  radar  data.  It  is 
designed  as  a  real  function  in  FORTRAN  that  accepts  time 


as  the  input  parameter  and  returns  the  exponentially  in¬ 
terpolated  altitude  from  the  radar  data.  The  radar  data 
are  taken  from  an  ASCII  file,  thus  allowing  easy  correc¬ 
tion  of  the  data. 

7.3.5  OTHER  BOOK-KEEPING  ROUTINES 

Other  routines  used  in  Pass  3  provide  the  follow¬ 
ing  functions: 

-  inserting  calculated  results  into  the  header; 

-  printing  out  results  on  the  console; 

-  printing  out  results  into  an  ASCII  file  for  quick 
reference; 

-  reseting  all  parameters  to  their  initial  values. 

7 . 4  UTILITIES 

A  number  of  utility  programs  have  either  been 
written  or  adapted  to  provide  additional  services  to  the 
user  of  the  reduction  package.  These  routines  include: 

plotting  data  on  the  Tektronix  603  storage  scope  for 
quick  viewing  after  reduction  has  been  made; 

extracting  conductivity  and  altitude  values  from  the 
header  of  an  already  reduced  flight  into  Basic-Plus 
virtual  array  (disk  file)  format  for  plotting  on  the 
PDP  11/45; 

printing  out  information  from  the  headers,  including 
launch  date,  launch  site,  physical  parameters  of  the 
probe,  time,  altitude  and  conductivities; 

transfering  files  (including  programs  and  data) 
between  RT  —  11  and  UNIX  with  the  program  "rtpip"; 
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creating  data  libraries  and  printing  out 
with  the  program  "PIP"; 

using  other  system  programs  from  RT-11  or 


programs 

UNIX. 


CHAPTER  VIII 


MICRO-PROCESSOR  BASED  SYSTEM 

The  data  reduction  system  is  designed  with 
micro-processor  adaptability  in  mind.  In  this  chapter, 
two  possible  micro-processor  implementations  are  dis¬ 
cussed.  One  is  at  the  system  level  using  Original  Equip¬ 
ment  Manufacturer  (OEM)  products  while  the  other  is 
designed  at  the  discrete  component  level. 

8 SYSTEM  APPROACH 

This  proposed  micro-computer  system  is  based  on 
the  PDP  11/03  micro-computer,  which  is  built  around  a 
LSI-1  1™  micro-computer  [Microcomputer  Handbook  (1976)]. 
The  LSI-11  is  the  smallest  of  the  PDP  11  family,  but  it 
shares  almost  the  same  architecture  and  instruction  set 
as  the  rest  of  the  PDP  11  processors  (Figure  8.1).  The 
LSI-11  interfaces  to  the  other  peripherals  through  the 
Q-BUSTM,  which  is  very  similar  to  the  UNIBUS.  The  inter¬ 
rupt  structure  on  the  LSI-11  is  the  same  as  on  the  PDP 
11,  and  therefore  the  same  data  acquisition  method  can  be 
used.  (However,  on  the  LSI-11  only  one  level  of  priority 
interrupt  is  present  instead  of  seven  as  on  the  PDP  11's. 
Our  system  only  requires  one  level.)  The  block  diagram  of 
this  system  is  shown  in  Figure  3.2  and  the  respective 
functions  will  now  be  discussed. 
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Using  ihe  now  11/03-LK  9-slot  back- 
p'une  version  of  The  POP  11/03  Ihe 
PDP  11T03-L  System  features  the 
new  RLV1 1  Disk  Subsystem  with  tv.o 
RlOI  5  2  mitt, on  byte  cartridge  disk 
dr. v:'S  tO-4  million  bytes) 

CONFIGURATION  DATA: 

Hardware  The  system  includes  the 
POP  11/03-LK  microcomputer  with 
32K  words  (6 A K  bytes)  RAM.  KEV11 
Extended  Arithmetic  Option  BDV11 
Eoctstrjp.  two  RLOI  Disk  Drives  and 
Confoller,  H9512  Cabinet,  and  your 
choice  ot  LA36  DECwnter  It  VT100 
CRT  Terminal,  or  LS120  DECwriter  III 
console  device  with  EIA  serial  inter¬ 
lace 

Software:  RT-l  1  Real-Time  Operat¬ 
ing  System  and  a  complete  set  of 
diagnostic  programs  for  this  configu¬ 
ration 


1  PDP  11/03  System  Configuration 


The  PDP  11/03  system  is  configured  with  32K  words 
of  MOS  memory.  The  processor  also  has  four  DLV11-J  asyn¬ 
chronous  serial  interfaces  to  provide  RS-232C  interface 
to  four  peripheral  devices,  including  the  LS-120  console 
terminal.  In  the  proposed  system,  a  Hewlett-Packard 
7221A  plotter  is  included  to  provide  hard  copy  outputs. 
A  dual  RL01  hard  disk  system  or  dual  RXV-21  double  densi¬ 
ty  floppy  disk  system  is  used  for  mass  storage  of  data 
and  programs.  The  RL01  hard  disk  can  provide  10  million 
bytes  of  on-line  storage  and  a  much  faster  data  transfer 
rate.  The  R XV 11—21  floppy  disk  system,  on  the  other 
hand,  can  only  provide  1  million  bytes  of  on-line 
storage,  but  the  low  cost  of  floppy  disks  makes  it  an  at¬ 
tractive  alternative. 

The  LPS  is  not  supported  by  the  LSI-11  processor 
but  a  number  of  modules  are  available  which  can  perform 
the  functions  of  the  LPS.  They  are  listed  below: 

-  KWV11  programmable  real  time  clock; 

-  ADV11  16  channel,  12  bit  A/D  converter; 

-  AAV11  4  channel,  12  bit  D/A  converter; 

-  DRV11  Digital  I/O  interface. 

The  hardware  descriptions  for  these  modules  and  other 
necessary  peripherals  are  listed  in  Table  8.1. 

The  Tektronix  603  storage  scope  is  used  for 
viewing  and  scaling  the  waveforms.  The  AAV11  D/A  con¬ 
verter  drives  the  storage  scope.  Both  the  storage  scope 


61 


FDII-.'IA  !  ''.'-11/2  Central  Processor  Unit  with  power  fp.il/ aulo  restart  LSI-11  bus  interface  and  vector  in- 

.....  n..r.c.inuj.  Sin.jie  board,  size:  8.9  in.  x  5.2  in.  (22.8cm  x  13.2cm). 

X.C .  11-DC  32K  byte  it, in  Com  Access  Memory  (RAM)  including  on-board  refresh.  Single  board,  size:  8.9  in.  x 

5.2  in.  (22.3cm  x  13.2cm). 

RLV11-AK  5  2  1.1  byte,  top  loading  removable  cartridge  subsystem  for  all  11/03L  LSI-11  systems,  two  quad 

control  boards,  size:  8  9  in  x  10  in.  (22  8cm  x  25  4cm)  NOTE:  Requires  H9273  or  BA11N  backplane. 

PL01-AK  Add  on  5  2  M  byte  top  loading  removable  cartridge  disk  drive. 

RXV21-L.  \  Dual  drive  dr  idle  density.  102-lK-bytos  capacity  plus  DMA  LSI-11  interface  board.  Single  bccru. 

RX  /21-3D  size:  8.9  in.  x  5.2  in.  (22.8cm  x  13  2cm). 

DLV11-J  F  ur  nc'upendcntly  configurable  serial  line  units.  Supports  RS-422  and  RS-423  (compat.ble 

T  -c32  C).  Selectable  parity,  data  and  stop  bits.  Baud  rates  Irom  150  to  33400.  Single  boarc.  „.ze: 

£...  in.  x  5.2  in.  (22.8cm  x  13.2cm).  NOTE:  Requires  one  cable  per  line  and  H3270  test  connectors  for 
c. -gnostics. 

LA120-BA  DECwriter  III  keyboard  printer  with  numeric  pad  which  operates  on  RS-232C  standard.  132-coI- 
l  in.  7  x  7  dot  matrix,  130-cps  smart  bi-directional  printing  designed  for  serial  1200  baud  com- 
:....r.ications.  20mA  current  loop  interlace  optional. 

1 -A  •1-chr.rncl  (buffeted),  D'A  converter  Output  ranges:  ±  2.56V,  ±5.12V,  ±10.24V,  C-5.12V, 

Id. 24V.  Single  board,  size:  8.9  in.  x  10  in.  (22.8cm  x  25.4cm). 

.'  OV11-A  12-1'  f.  IB-rhanncl.  single-ended.  (6  ch  .ar.cl  differential)  A/D  converter.  Input  ranges:  ^5  12V. 

.  ..  circuitry.  50/is  conversion  time.  Single  board,  size:  8.9  in.  x  10  in.  (22.0em  x 

'.)■ 

KVVV11-A  Programmable  Crystal  Cloc'^with  frequencies  from  100  Hz  to  1  MHz  plus  60  cycle  and  external 

ir .it.  Single  board,  size:  8.9  in.  x  10  in.  (22.8cm  x  25.4cm). 

DF-.Y'il  P  'cl  Line  Interface  Unit.  16-bit  diode-damped  input;  16-bit  latched-drive  output.  Protocol  and 

...-cl  signals.  Suggested  cables  BC07D  or  BC03R.  BC08R  can  be  used  as  a  diagnostic  cable. 
Single  board,  size:  6.9  in.  x  5.2  in.  (22.8cm  x  13  2cm). 

BA- "  ;  E'-.r.a.nr’or  Box.  Includes  4x9  b.-  ikplano  (LSI-11  bus  on  slots  A  &  B  only)  with  240  WATT  power 

LA,  i -13-  sir  ,,!y. 

P.T-11  F.  '-Vino  Operating  System  designed  for  the  single  user.  Includes  single  |Ob  monitor  and  foreground/back¬ 
ground  m.  tor  RT-11  system  programs  include  EDIT.  MACRO-11,  LINKER.  ODT,  and  utilities.  Minimum  hardware: 

CPU  witn  16K  by.o  (32K  byte  for  foreground/background  operation). 

Par'..  -11  '  m  ien  P.l.ind.wd  Fortran  with  symbolic  debugger  usable  on  LSI  Bus  products,  also  includes  a  li¬ 
brary  oi  .  ..ere,.  .  :■  b/oulmcs  not  usable  on  LSI  Bus  products.  Prerequisite:  license  to  use  RT-11. 


Table  8.1  LSI -11  System  Supporting  Hardware  Description 


and  the  D/A  converter  have  12-bit  resolution  (4096)  so 
they  interface  nicely.  The  Hewlett-Packard  plotter,  with 
the  proper  software  configuration,  can  draw  waveform  im¬ 
ages  on  to  a  hard  copy  in  addition  to  plotting  the  re¬ 
duced  conductivity  data.  This  would  permit  the  detailed 
study  of  waveforms  and  probe  response  characteristics  at 
a  higher  resolution  than  what  is  presently  possible.  The 
Schmitt-trigger  available  on  the  LPS  can  very  easily  be 
built  using  op  amps  to  trigger  the  KWV11  programmable 
real-time  clock.  If  a  second  KWV11  is  available,  it  can 
be  used  in  place  of  the  DRV11  digital  I/O  and  the  time 
code  generator.  It  can  be  programmed  to  run  at  any 
desired  frequency  for  timing  information.  For  example, 
10  Hz  will  give  0.1  second  resolution;  100  Hz  will  give 
0.01  second  resolution  in  timing.  In  order  to  prevent 
overflow  from  the  16  bit  word ,  the  time  count  can  be 
stored  in  2  words  using  the  "ADD"  and  "ADC"  instructions. 
In  this  case,  for  example,  up  to  42949673  seconds  can  be 
stored  with  0.01  second  accuracy'!’'. 

The  PDP  11/03  system  has  the  RT-11  V3B  operating 
system  and  FORTRAN  IV  software  support.  Plotting 
subroutines  written  in  FORTRAN  are  also  available  from 
Hewlett-Packard.  Therefore,  this  system  could  be  operat¬ 
ing  in  a  very  short  period  of  time  with  minimum  software 

''  If 'the  clock  frequency  =  100  Hz,  i.e.  each  second  has 
100  counts,  then  32  bits  provides  2 32 _ i  counts;  there¬ 
fore  time  =  (23 2  _  1  ) / 1 00  =  42949673  seconds. 


modifications  to  our  existing  system. 

8.2  COMPONENT  APPROACH 

Due  to  the  enormous  number  of  micro-processor 
chips  available,  this  discussion  will  not  attempt  to 
focus  on  any  particular  chip.  Instead,  a  general  confi¬ 
guration  will  be  outlined  so  that  the  system  designer  can 
use  it  as  a  guide  line  for  the  particular  design  (Figure 
8.3). 

The  heart  of  the  system  is  the  data  acquisition 
phase  as  described  in  this  paragraph.  In  the  component 
approach,  a  fixed  frequency  clock  is  fed  into  a  binary 
counter  for  frequency  measurement.  When  a  pulse  is 
detected  from  the  incoming  signal,  the  Schmitt- trigger  is 
fired  causing  the  contents  of  the  counter  to  be  latched 
into  the  tri-state  latch.  It  also  sets  the  flip-flop 
that  is  tied  to  the  interrupt  request  line  (INTR)  of  the 
micro- processor  chip.  (A  flip-flop  is  used  because  most 
micro-processors  will  only  sample  the  interrupt  request 
line  after  the  last  machine  cycle  of  each  instruction. 
If  the  flip-flop  were  not  used,  the  interrupt  request 
might  be  lost  [Liu  (1979)3.)  After  the  (INTR)  signal  is 
recognized  by  the  micro-processor,  an  interrupt  ack¬ 
nowledge  (INTA)  signal  will  be  issued  to  reset  the  flip- 
flop.  Depending  on  the  hardware  and/or  software  confi¬ 
guration  of  the  particular  micro-processor  system,  con- 


Figure  8.3  Micro-processor-based  Data  Acquisition  System 
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trol  is  passed  in  someway  to  an  interrupt  service 
routine.  The  interrupt  service  routine  will  issue  an  I/O 
read  to  read  the  content  of  the  counter.  This,  in  turn, 
causes  the  read  (RD)  line  to  activate  the  output  enable 
(OE)  line  on  the  tri-state  latch  which  puts  the  count  on 
the  data  bus  of  the  micro-processor  system.  The  count 
can  therefore  be  read  in  from  an  I/O  port  of  the  micro¬ 
processor.  (The  address  of  the  port  is  determined  by  the 
I/O  decoding  logic.)  After  the  interrupt  service  routine 
finishes  its  task,  control  returns  to  the  normal  program 
sequence  and  a  ready  (RDY)  signal  will  be  issued.  This 
signal  resets  the  counter  which  is  ready  for  the  next 
pulse. 

After  the  data  are  digitized,  they  need  to  be 
stored  and  reduced  using  either  a  micro-processor 
development  system  or  a  mini-computer. 

8.3  COMPARISON  BETWEEN  THE  TWO  APPROACHES 

The  system  approach  is,  of  course,  more  costly. 
Yet  the  reliability  and  maintenance  available  from  OEM 
products  impose  a  much  greater  advantage  over  the  com¬ 
ponent  approach.  In  addition,  the  compatibility  of  the 
LSI-11  and  PDP  11  computer  systems  allows  minimum  modifi¬ 
cations  to  the  existing  software. 

The  component  approach  requires  a  tremendous 
amount  of  design  effort  but  with  less  hardware  cost. 


Therefore,  it  can  also  be  an  attractive  alternative  for 
production  environment. 


CHAPTER  IX 


ERROR  ANALYSIS 

In  this  chapter,  errors  that  might  possibly 
develop  during  the  different  stages  in  the  data  reduction 
process  are  considered  and  their  significances  to  the 
reduction  results  are  analyzed.  Improvements  to  the 
present  system  are  proposed. 

9-1  ERRORS  RESULTING  FROM  A  FIXED  RATE  CLOCK  AND  A  FIXED 
WORD  LENGTH 

Two  kinds  of  errors  result  from  using  a  fixed 
frequency  clock  reference  on  a  fixed  word  length  comput¬ 
er.  The  first  one  is  very  obvious.  Since  the  clock  is 
counting  at  100  kHz,  any  incoming  signal  with  a  frequency 
greater  than  100  kHz  is  not  detected.  The  normal  operat¬ 
ing  frequency  of  the  probe  is  between  0  and  200  Hz  only 
and  thus  the  significance  of  this  kind  of  error  is  very 
small. 


The  other  kind  of  error  results  from  the  digital 
computer  having  a  fixed  word  length.  In  the  case  of  the 
PDP  11 's,  the  word  length  is  16  bits.  This  means  that  an 
unsigned  integer  representation  can  at  most  be  2^°-1, 
i.e.  65535.  In  other  words,  when  the  clock  is  counting 
at  100  kHz,  only  65535  counts  can  be  recorded  before 
overflow  occurs.  This  number  corresponds  to  0.65535 
second  or  a  frequency  of  1.53  Hz.  In  normal  operation,  a 
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frequency  lower  than  1.53  Hz  is  unlikely  but  this  situa¬ 
tion  is  possible  if  there  is  signal  dropout.  In  this 
case,  the  clock  counts  up  to  65535  counts,  and  then  any 
count  beyond  that  will  reset  the  counter  to  zero  and  it 
starts  counting  from  zero  again.  In  effect,  0.65535 
second  is  lost  and  the  count  at  the  end  of  this  pulse  is 
not  correct.  This  was  a  major  problem  in  the  past  and 
there  was  no  easy  solution  because  of  the  chosen  design 
model.  In  the  present  system,  however,  the  problem  is 
insignificant  because  the  introduction  of  the  external 
time  code  generator  has  completely  separated  the  pro¬ 
grammable  real-time  clock  from  the  task  of  keeping  the 
in-flight  timing.  Whenever  dropout  occurs,  only  one  bad 
point  is  introduced.  The  timing  distortion  associated 
with  this  point  is  not  carried  beyond  that  particular 
waveform  because  each  time  the  push  button  is  pressed, 
Pass  1  will  go  out  to  the  external  time  code  generator  to 
get  timing  information  which  re-synchronizes  itself. 

When  signal  dropouts  occur  during  the  data  ac¬ 
quisition  phase,  the  operator  can  allow  the  dropouts  to 
elapse  until  normal  signal  playback  recovers .  Then  the 
push  button  switch  can  be  issued  and  timing  is  re¬ 
synchronized.  During  the  data  reduction  pass,  the  opera¬ 
tor  will  notice  that  the  waveform  is  shortened  due  to  the 
data  loss  during  dropout.  The  waveform  can  be  discarded 
and  reduction  can  resume  on  the  next  waveform. 
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9.2  ERRORS  RESULTING  FROM  I/O  ACTIVITIES 

Pass  1  is  the  time-critical  phase  since  it  has  to 
capture  the  incoming  signal  in  real  time.  Looking  at  the 
data  acquisition  method  discussed  before,  the  time- 
critical  portion  of  Pass  1  is  the  interrupt  service 
routine.  When  the  Schmitt- trigger  is  fired  by  an  incom¬ 
ing  pulse,  an  interrupt  is  generated.  Control  is  then 
passed  to  the  interrupt  service  routine  and  during  this 
period,  further  interrupts  are  disabled.  If  another 
pulse  comes  in  during  this  time,  it  is  not  recognized. 
Coding  of  che  interrupt  service  routine  has  been  careful¬ 
ly  optimized  to  reduce  execution  time  (see  Appendix  D.1). 
With  our  particular  hardware  and  software  configuration, 
under  normal  circumstances,  each  pulse  takes  125.5  us!2 
to  be  stored.  The  corresponding  bandwidth  is  7.95  kHz 
(including  the  A/D  converter  sampling  time).  In  the 
worst  case  for  which  the  buffer  has  been  filled  and  I/O 
activity  is  required,  215.9  us  are  needed,  thus  giving  a 
bandwidth  of  6.35  kHz.  With  the  normal  operating  range 
of  the  probe  frequencies  (0  to  200  Hz) ,  this  error  is  in¬ 
significant. 

The  above  calculations  are  based  on  the  fact  that 
the  RT-11  operating  system  has  the  full  support  of  asyn¬ 
chronous  I/O.  Each  time  I/O  activity  is  needed,  the  pro¬ 
gram  only  issues  the  I/O  request  through  the  use  of  "sys- 
^  1  us  =  10~®  second. 
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tem  macros”.  Control  then  returns  to  the  calling  program 
after  the  I/O  request  has  been  recognized.  The  I/O  ac¬ 
tivities  will  be  executed  in  the  background  asynchronous¬ 
ly  [RT-11  System  Reference  Manual  (1975)]. 

9.3  ERRORS  RESULTING  FROM  DISPLAY  AND  CURVE  FITTING 

In  the  field  of  computer  graphics,  there  has  al¬ 
ways  been  the  problem  of  mapping  the  world  coordinates 
into  the  device  coordinates  since  the  world  coordinate 
system  has  almost  unlimited  addressable  space  whereas  the 
device  coordinate  system  is  always  restricted  by  the  lim¬ 
ited  discrete  representation  of  physical  addressable 
spaces.  The  same  problem  arises  here  since  the  count 
from  the  programmable  real-time  clock  can  range  from  1  to 
65535,  i.e.  there  are  65536  possible  readings.  However, 
there  are  only  4096  addressable  points  on  the  storage 
scope,  thus  implying  that  points  with  close  values  will 
fall  into  the  same  spot  on  the  screen.  Since  the  resolu¬ 
tion  of  the  human  eye  is  far  less  than  the  scope's  reso¬ 
lution,  this  effect  is  not  significant. 

The  next  problem  arises  from  the  basic  relation¬ 
ship  of  calculating  frequency  from  count.  These  two 
parameters  are  inversely  proportional  to  each  other  and 
thus  a  hyperbolic  relation  is  expected.  In  contrast  to  a 
linear  relationship  between  count  and  frequency,  the  fre¬ 
quencies  in  the  middle  portion,  for  example  from  50  to 
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160  Hz,  will  only  be  resolved  by  1,041  counts^3.  This 
means  a  1-in-10  resolution  instead  of  1-in-100,000  as  ex¬ 
pected.  This  in  part  is  the  reason  for  choosing  such  a 
high  reference  frequency  as  100  kHz. 

The  other  type  of  error  comes  from  the  least 
square  curve  fitting  program.  The  algorithm  for  slope 
fitting  is  designed  such  that  a  weighting  factor  between 
0.0  and  1.0  is  assigned  to  each  point  according  to  how 
far  that  point  is  from  the  first  fitted  line.  The  impli¬ 
cation  is  that  points  farther  away  from  the  fitted  line 
will  bear  less  significance  to  the  final  fit.  However, 
there  are  two  important  points  to  remember.  First,  the 
choice  of  the  time  segment  to  be  scaled,  i.e.  the  two  end 
points,  is  very  important.  If  the  end  points  are  not 
correct,  a  correct  result  is  less  likely  to  occur. 
Secondly,  if  there  exists  an  overwhelming  amount  of  noisy 
data  in  the  region  between  the  two  end  points,  the  resul¬ 
tant  slope  will  not  be  correct.  This  is  largely  because 
the  first  trial  will  assign  an  equal  weight  to  every 
data  point.  If  the  first  trial  is  not  correct,  as  in  the 
case  of  a  lot  of  noise,  the  successive  trials  will  not 
possioly  be  correct.  To  solve  this  problem,  a  small 
clean  portion  of  the  waveform  needs  to  be  chosen  to  gen¬ 
erate  the  fit.  This  usually  achieves  a  very  satisfactory 
result. 

TTIT _cfvtl  ~ TUiJkHz/60  =~TC>66;  cnt2  =  100kHz/ 160  =  625; 
then  (cntl  -  cnt2)  =  1000. 
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Finally,  the  reduction  of  waveforms  having  very 
steep  slopes  (large  conductivity  values)  requires  the  ex¬ 
pansion  routine  to  blow  up  a  very  small  portion  of  the 
waveform  on  the  screen.  This  allows  pin-pointing  to  the 
exact  end  points  and  in  turn,  scaling  of  the  proper  time 
segment  of  the  waveform. 

9.M  ERROR  RESULTING  FROM  THE  TIME  CODE  GENERATOR 

For  the  present  hardware  available,  the  time  code 
generator  is  clocked  every  second.  Therefore,  timing  can 
be  at  most  one  second  off  from  the  actual  value.  There 
are  two  ways  to  improve  this  accuracy.  The  first  method 
requires  a  time  code  generator  that  gives  more  resolution 
in  output  timing.  The  second  method  is  to  put  in  a 
correction  factor  to  include  the  timing  elapsed  up  to  the 
interval  of  the  designated  slope.  This  can  very  simply 
be  done  by  summing  up  all  the  counts  to  the  middle  point 
of  the  time  interval  scaled.  By  storing  this  corrected 
time  factor  in  the  extra  space  of  the  header,  the  timing 
accuracy  is  improved.  This  would  be  especially  important 
during  the  early  stage  of  the  flight  where  the  payload  is 
decending  at  a  relatively  high  velocity. 


CHAPTER  X 


RESULTS  AND  FUTURE  IMPROVEMENTS 

In  this  chapter,  electrical  conductivity  values 
obtained  from  the  computerized  reduction  system  are 
presented.  In  addition,  a  section  is  presented  describ¬ 
ing  some  of  the  problems  possibly  encountered  during  the 
reduction  process  and  how  they  might  be  corrected.  Fi¬ 
nally,  possible  hardware  and  software  improvements  are 
discussed  . 

K).J_  RESULTS 

The  end  results  of  this  computerized  reduction 
scheme  are  the  electrical  conductivity  values  measured  in 
the  middle  atmosphere.  Using  this  reduction  procedure, 
conductivity  data  from  four  blunt  probes  recently 
launched  on  super  Loki  rockets  are  shown  in  Figures  10.1 
to  10.4.  All  four  rocket  launches  were  conducted  at  re¬ 
latively  high  latitudes.  The  two  1977  rocket  launches 
occurred  at  Poker  Flat,  Alaska  (65°  N,  147°  W)  as  part  of 
a  program  to  study  the  ionization  effects  of  auroral  en¬ 
ergetics  on  the  middle  atmosphere.  The  two  1979  launches 
were  part  of  a  solar  eclipse  rocket  program  conducted  at 
Red  Lake,  Ontario,  Canada  (51°  N,  94°  W) . 

In  all  of  the  figures,  the  plus  and  minus  signs 
represent  positive  and  negative  conductivity  values, 
respectively,  while  the  dots  indicate  where  the  measure- 
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ments  are  comparable  in  value.  The  daytime  measurements 
(Figure  10.1,  10.2  and  10.4)  at  higher  altitudes  show  no¬ 
ticed  differences  in  the  polar  conductivity  components 
for  the  same  altitude,  with  the  negative  conductivity 
values  being  relatively  larger.  The  larger  negative  con¬ 
ductivity  values  are  thought  to  demonstrate  the  presence 
of  electrons,  which  have  relatively  larger  mobility 
values.  In  general,  more  variability  is  observed  in  the 
reduced  negative  conductivity  measurements,  which  is  in¬ 
dicative  of  the  difficulty  in  reducing  the  electron 
current  component. 

The  nighttime  conductivity  measurements  (Figure 
10.3)  generally  show  little  differences  in  the  polar  con¬ 
ductivity  components  at  the  same  altitude,  which  is  ex¬ 
pected  when  there  are  no  electrons  present.  Surprising¬ 
ly,  the  negative  conductivity  values  at  lower  altitudes 
are  relatively  smaller  than  the  corresponding  positive 
values,  which  is  generally  not  the  observed  case.  These 
relatively  smaller  negative  conductivity  values  are  prob¬ 
ably  attributed  to  smaller  ion  mobility  values. 

U).2  USER'S  EXPERIENCE 

The  advantages  of  this  data  reduction  package  are 
enhancement  in  speed,  greater  flexibility  and  improved 


accuracy . 
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CONDUCTIVITY  Cnho/cn) 

Electrical  Conductivity  Measurements  for  February  26,  1979  at  2240  CST 


Electrical  Conductivity  Measurements  for  February  27,  1979  at  0840  CST 
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Digitization  only  needs  to  be  done  once  for  each 
flight  and  therefore  a  lot  of  play-back  time  is  saved. 
This  system  frees  the  user  from  tedious  manual  scaling  of 
waveforms.  Scaling  (by  a  numerical  method),  line  genera¬ 
tion,  calculation  and  data  storage  are  all  done  automati¬ 
cally.  Furthermore,  the  results  are  all  self-contained 
in  the  header  for  each  waveform  of  each  flight.  This  in¬ 
formation  is  in  machine  readable  form  and  can  be  utilized 
directly  by  machines.  For  example,  a  plotting  utility 
program  has  been  written  to  read  the  results  from  the 
headers  and  plot  them  out  on  the  Tektronix  603  storage 
scope  for  viewing.  Another  routine  has  been  written  to 
extract  data  from  the  header  and  plot  the  results  on  the 
hard  copy  plotter  of  the  PDP  11/45.  Different  outputs 
are  thus  generated  according  to  the  format  required 
without  any  complications.  Different  sizes  of  plots, 
different  titling  and  labels  are  all  done  easily. 

Each  45-minute  flight  occupies  approxmately  1000 
blocks  of  disk  space.  If  storage  is  tight,  the  Schmitt- 
trigger  level  on  the  incoming  signal  can  be  turned  off 
periodically  during  the  latter  part  of  the  flight  where 
the  altitude  changes  slowly.  Since  the  external  timer 
runs  asynchronously,  timing  is  still  kept  correctly. 

Accuracy  is  another  reason  for  using  the  reduc¬ 
tion  package.  The  data  digitization  method  enhances  the 
resolution  and  in  addition,  the  expansion  routines  can 
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blow  up  any  portion  of  the  waveform  so  that  detailed 
study  of  the  waveform's  characteristics  are  possible. 
Using  numerical  techniques  to  compute  the  slope  are  also 
advantageous,  especially  for  very  steep  or  very  flat 
slopes . 

Finally,  sporadic  noise  pulses  or  telemetry 
drop-outs  will  often  corrupt  data  waveforms  displayed  on 
a  strip  chart  such  that  they  are  not  suitable  for  scal¬ 
ing.  It  has  been  observed  that  these  effects  are  often 
reduced  using  computer  techniques  such  that  conductivity 
information  can  still  be  extracted  from  the  waveform. 

J_0 . 3  FUTURE  IMPROVEMENTS 

In  this  section,  future  improvements  to  the  data 
reduction  system  are  discussed.  The  improvements  include 
steps  toward  a  more  automated  reduction  system,  more  so¬ 
phisticated  hardware  and  more  fully  coordinated  software. 

10.3*1  TOWARDS  AUTOMATION 

Probably  the  most  significant  step  toward  automa¬ 
tion  of  the  reduction  process  actually  involves  a  probe 
redesign  such  that  an  indication  of  the  time  for  zero- 
volt  probe  potential  is  also  telemetered  with  the 
cur  rent-voltage  response.  Such  an  indication  would  per¬ 
mit  adaptation  of  a  new  computerized  waveform  segmenta¬ 
tion  scheme  to  replace  the  present  operator-controlled 
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procedure.  In  addition,  it  would  help  in  identifying  the 
two  separate  polar  currents  collected  by  the  probe.  Fi¬ 
nally,  it  could  be  used  to  estimate  probe  potential  dur¬ 
ing  other  times  of  the  waveform,  which  is  particually 
relevant  in  the  further  reduction  of  probe  data  to  esti¬ 
mate  electron  density  (for  blunt  probes)  and  ion  density 
(for  Gerdien  condensers). 

20.3.2  HARDWARE  IMPROVEMENTS 

More  sophisticated  hardware  will  of  course,  make 
the  data  reduction  system  more  attractive.  A  display 
processor  equipped  with  a  light  pen  could  replace  the 
storage  scope,  thus  allowing  a  more  convenient  way  of 
scaling  waveforms.  A  plotter  interfaced  to  the  system 
not  only  would  provide  a  final  form  of  output  but  with  a 
proper  software  setup,  would  draw  each  scaled  waveform  on 
paper  as  well.  This  would  provide  an  opportunity  for 
studying  the  detailed  structure  of  the  probe's  current 
response  characteristics.  Lastly,  a  second  programmable 
real-time  clock  could  be  used  for  obtaining  accurate  tim¬ 
ing  information  under  software  control. 

2JD.3.3  SOFTWARE  IMPROVEMENTS 

With  the  digitized  data  entered  into  the  comput¬ 
er,  future  software  improvements  are  concerned  with 
further  manipulation  of  this  stored  information.  Specif¬ 
ically,  further  software  development  might  include:  im- 
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provement  of  the  method  for  accumulating  time;  estimation 
of  probe  potential;  calculation  of  electron  number  densi¬ 
ty;  refinement  of  plotting  routines  to  include  multi¬ 
images  on  each  output;  and  expansion  of  library  routines. 
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APPENDIX  A 


PROGRAMS  USED  IN  THE  DATA  REDUCTION  SYSTEM 
A.J_  PROGRAMS  UNDER  THE  RT-JJ_  SYSTEM 
PASS1  -  Data  acquisition 

PASS2  -  Data  segmentation  and  header  insertion 
PASS3  -  Data  display  and  reduction 
PLOT  -  Plotting  conductivity  values  vs  altitude  on 
the  Tektronix  603  storage  scope 
XTRACT  -  Extract  conductivity  values  from  the  header 
and  put  them  into  Basic-Plus  virtual  array 
format  for  plotting  on  the  PDP  11/45 
PIP  -  DEC'S  system  program  for  printing  or 

transfering  files  between  peripheral  devices 

A. 2  PROGRAMS  UNDER  THE  UNIX  SYSTEM 

rtpip  -  Transfering  files  between  the  UNIX  and  RT-11 
systems 

plot  -  General  purpose  plotting  routine 
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APPENDIX  B 


USER'S  MANUAL 


B._l  HARDWARE  SETUP  FOR  PASS 

[i]  See  the  schematic  in  Figure  4.1. 

[ii]  Connect  the  output  of  the  playback  tape  recorder  to 
the  input  of  the  Schmitt- trigger  2  on  the  LPS. 

[iii]  Connect  the  pulse  generator  to  channel  0  of  the  A/D 
converter  on  the  LPS. 

[iv]  Connect  the  terminating  signal  output  to  channel  2 
of  the  A/D  converter  on  the  LPS. 

[v]  If  there  is  another  data  signal,  connect  it  to  chan¬ 
nel  1  of  the  A/D  converter. 

[vi]  A  strip  chart  machine  can  optionally  be  connected  to 
the  output  of  the  Schmitt-tr igger  2  to  obtain  a  hard 
copy  of  the  data  waveform. 

[vii]  An  oscilloscope  can  optionally  be  connected  to  both 
the  input  and  the  output  of  the  Schmitt- trigger  2 
for  monitoring. 

B . 2  OPERATING  PROCEDURE  FOR  PASS  1 

[i]  Set  up  the  hardware  configuration  as  described  in 
Append ix  B. 1 . 

[ii]  Turn  on  all  the  equipment. 

[iii]  Play  back  portions  of  the  tape  for  a  few  times  to 
determine  the  correct  triggering  level  on  the 
Schmitt-tr igger  2. 

[iv]  Reset  the  time  code  generator. 

[v]  Boot  up  RT  —  11  and  execute  the  program  "PASS1.SAV". 

[vi]  Erase  the  storage  scope. 

[vii]  Start  the  tape  recorder.  The  storage  scope  will  in¬ 
dicate  the  data  are  entered. 

[viii]  Start  the  time  code  generator  when  launch  is 
detected  (either  from  the  scope  or  the  strip  chart). 
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[ix]  Push  the  switch  on  the  pulse  generator  whenever  a 
new  waveform  is  recognized. 

[x]  Issue  the  terminating  signal  (a  5  V  falling  pulse) 
when  the  data  end. 

[xi]  Execute  the  program  "PIP"  to  see  if  the  file 
"PASS  1 . TMP"  is  present. 

B.3  OPERATING  PROCEDURE  IN  PASS  2 

[i]  Obtain  data  from  Pass  1  as  described  in  Appendix 
B.2. 

tii]  Boot  up  RT-11 . 

[iii]  Execute  the  program  "PASS2.SAV". 

[iv]  Answer  all  the  questions  the  program  asks  concerning 
the  probe's  information. 

[v]  The  program  will  execute  automatically.  The  LED 

display  on  the  LPS  will  show  the  time  of  each 

waveform  sequentially. 

[vi]  When  the  program  terminates,  the  file  "PASS1.TMP" 
can  be  deleted  using  "PIP". 

[ vii ]  Use  "PIP"  to  check  if  "PASS2.TMP"  is  present. 

B. 4  OPERATING  PROCEDURE  IN  PASS  3 

[i]  Turn  on  the  storage  scope  and  erase  it. 

tii]  Boot  up  RT-11  and  execute  "PIP"  to  find  out  the  size 

of  the  file  "PASS2.TMP". 

t iii ]  Execute  "PASS3 • SAV" . 

tiv]  Answer  the  questions  requested  by  the  program: 

a)  The  name  of  the  altitude  file  is  an  ASCII  file 
that  contains  the  time  vs  altitude  radar  data. 

b)  The  number  of  blocks  in  "PASS2.TMP"  is  the  size 
obtained  in  step  [ii]  above. 

c)  The  calculation  constant  is: 
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FOR  BLUNT  PROBE 


for  conductivity  calculations. 

[v]  The  program  proceeds  and  asks  from  which  waveform 
number  to  begin  reduction.  Reduction  can  start  at 
any  waveform  desired.  A  <CR>  implies  the  first. 

[vi]  The  complete  waveform  chosen  will  appear  on  the 
screen  with  a  vertical  scale  from  0  to  200  Hz.  The 
time  for  that  waveform  will  appear  on  the  LED 
display  as  seconds  into  the  flight. 

[ v ii 3  The  program  will  ask  if  expansion  in  the  X  and/or  Y 
direction(s)  is(are)  desired.  If  no  expansion  is 
needed,  go  to  step  [ix]. 

[viii]  The  program  will  request  the  X  and  Y  boundaries 
for  expansion.  A  <CR>  alone  means  display  the  whole 
waveform  from  0  to  200  Hz  again.  The  boundary  sup¬ 
plied  in  the  X  direction  means  only  displaying  in 
the  selected  X-direction  limit  but  from  0  to  200  Hz 
vertically.  If  both  the  X  and  the  Y  directions  are 
given,  the  display  will  be  on  the  selected  portion 
in  both  directions.  The  program  will  then  go  back 
to  step  [vii]  and  repeat  until  no  more  expansions 
are  desired. 

[ix]  The  program  will  ask  for  the  two  end  points  on  the 
positive  conductivity  slope.  If  a  <CR>  is  given,  go 
to  step  Cxi]. 

[ x 3  The  program  receives  the  two  end  points  and  will 
find  the  best-fitted  straight  line  between  them. 
Then  a  solid  line  will  be  generated  on  the  screen. 
The  program  will  ask  for  break  points  again.  If  the 
fit  is  satisfactory,  type  a  <CR>,  otherwise  input 
another  two  points  and  step  [x]  will  repeat. 

[xi]  If  step  [x]  has  been  executed,  that  means  a  positive 
conductivity  value  is  present;  go  to  step  [xii], 
otherwise  go  to  step  [xiv]. 

[xii]  The  program  will  execute  step  [vii]  through  [xi] 
again  for  negative  conductivity  calculations.  If  a 
<CR>  is  given,  that  means  the  positive  and  negative 
conductivities  are  equal. 

[xiii]  The  program  will  ask  for  another  set  of  positive 
and  negative  conductivities.  If  they  are  present, 
steps  [vii]  through  [xiii]  will  repeat. 

[xiv]  The  reduced  conductivities,  altitude  and  time  will 
be  printed  on  the  console. 


[xv]  The  program  will  ask  to  continue  or  skip  forward.  A 
"C"  means  continue  to  the  next  waveform;  a  "Sw  means 
to  skip  some  of  them.  (The  program  will  ask  how 
many  waveforms  are  to  be  skipped.)  A  <CR>  means  to 
terminate  the  process. 

[xvi]  The  process  will  be  repeated  until  termination  is 
issued  from  step  [xv]  or  until  the  end  of  data  is 
encountered . 

[ xv ii ]  An  ASCII  file  called  "PASS3. TMPH  with  all  the  re¬ 
duced  information  has  been  created.  It  can  be  exam¬ 
ined  or  printed. 

B . 5  LIBRARY  CREATION 

[i]  A  data  library  can  be  created  by  using  PIP. 

[ii]  The  following  convention  for  extensions  on  data 
files  has  been  employed  to  distinguish  between  dif¬ 
ferent  file  contents: 

*.CAL  calibration  waveform  file  (PASS2  format) 
*.DAT  in-flight  data  file  (PASS2  format) 

* . LST  results  (ASCII) 

*.POS  positive  conductivity  vs 

altitude  file  (virtual  array) 

*.NEG  negative  (virtual  array) 

*.COM  positive  =  negative  (virtual  array) 

*.ALT  time  vs  altitude  radar  data  (ASCII) 


APPENDIX  C 


PROGRAM  LINKAGE 

The  following  commands  give  the  program  modules' 
names  and  linkage  procedure  for  producing  the  proper  exe¬ 
cutable  programs  under  the  RT-11  Operating  System.  They 
properly  compile  and  assemble  the  main  programs  and 
subroutine  3  into  object  modules  and  use  the  RT-11  linker 
(LINK)  to  produce  executable  programs  as  follows: 

PASS1 . SAV<PASS 1 . OB J 

PASS2.SAV<MAIN2.0BJ,PASS2. OBJ, QUEST. OB J/F 

PASS 3 . SAV<PASS3. OB J , SUB  1 . OBJ , SUB 2. OBJ , SUB 4 . OB J , SLOPE . OB J/F 

PLOT . SAV<PLOT . OB J/F 

XTRACT. SAV  <XTRACT. OB J , SUB  1 . OB J/F 


APPENDIX  D 


PROGRAM  LISTINGS 


u 
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PASS1  OF  DATA  PROCESSING 


THIS  ROUTINE  USES  THE  LPSKW  TO  COUNT 
AT  100  KHZ  AND  THE  SCHMITT -TRIGGER-2  TO 
COUNT  INCOMING  PULSE  FREQUENCY 
INPUTTED  FROM  THE  TAPE  RECORDER 


THE  ADCS  PROVIDE  COMMANDS  FOR  THE  ROUTINE 
AS  FOLLOWS: 

CHANNEL  0  —  A  +5V  PULSE  OF  50  MS  IS 
USED  FOR  WAVEFORM  SEGMENTATION 
TIMING  INFORMATION  WILL  BE  SAMPLED 
FROM  THE  DIO  PORTS  OF  THE  LPSDR 

CHANNEL  1  —  NOT  USED  AT  PRESENT  BUT  IS 
RESERVED  FOR  FUTURE  EXPANSION 

CHANNEL  2  —  TERMINATION  SIGNAL 
NORMAL  HIGH  (+5V) 


SOME  CONSTANT  DEFINITIONS 

TO  USE  HIGHLY  PARAMETERIZED  PROGRAMMING 
CAN  REDUCE  CHANCE  FOR  MAKING  ERROR  WHEN 
THE  PROGRAM  HAS  TO  ADAPT  TO  OTHER  SYSTEM 


BUFSIZ 

FILSIZ 

RECLEN 

CHANO 

CHAN1 

CHAN2 

KWENBL 

ERASE 

VCRDY 

EOWF  a 

NOINT 

ZEROV 

YVAL  = 

ICNT1 

ICNT2 

XINC  s 


=  4096. 
=  1000. 
=  16. 

1 

401 
1001 
=  1505 
:  10000 
:  2012 
0 

:  340 
:  5000 
4000 
=  4000. 

:  4000. 

2 


BUFFER  SIZE  =  4K 

DEFAULT  FILE  SIZE  =  1000  BLOCKS 

4K  BUFFER  =  16  BLOCKS  LONG 

BIT  PATTERN  TO  START  A  ADC  IN  CHO 

START  ADC  AT  CHANNEL  1 

START  ADC  AT  CHANNEL  2 

START  LPSKW.  100KHZ,  REPEAT  MODE 

ERASE  THE  SfcOPE 

LPSVC  Y-MODE.  FAST  INTENSIFY 

END  OF  WAVEFORM  MARK 

BR7,  NO  OTHER  INTERRUPT 

DIGITAL  VALUE  OF  0  V  FOR  ADC 

MIDDLE  OF  SCREEN 

COUNT  FOR  IDEL  LOOP1 

COUNT  FOR  IDEL  LOOP 2 

DISTANCE  BETWEEN  POINTS 


LPS  ADDRESSES  DEFINITION 

REFER  TO  LPS  OPERATING  MANUAL  FOR 

FURTHER  REFERENCES 

. NLIST 

ADST  a  170400  ;  LPSAD  STATUS 

ADBF  a  170402  ;  LPSAD  BUFFER 

KWST  a  170404  ;  LPSKW  STATUS 
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PASS1: 


KWBP  =  170406 
DR ST  =  170410 
DRIN  =  170412 
DROUT  =  170414 
VCST  =  170416 
VOX  =  170420 
VCY  =  170422 
KWIV  =  324 


LPSKW  BUFFER/PRESET 
DIGITAL  I/O  STATUS 
DIGITAL  I/O  INPUT 
DIGITAL  I/O  OUTPUT 
LPSVC  STATUS 
LPSVC  X 
LPSVC  Y 

LPSKW  INTERRUPT  VECTOR 
NOTE  THIS  IS  NON-STANDARD 


.MACRO 
JSR 
.  ENDM 

! MACRO 

RTS 

.ENDM 


MACRO 


fSTB 

BPL 

.ENDM 


SOME  MACRO  DEFINITIONS 

JUMP  TO  SUBROUTINE  THRU  PC 
CALL  A 
%7  ,A 

RETURN  FROM  SUBROUTINE 
RETURN 

n 

WAIT  IDLE  UNTIL  READY  BIT  (BIT  7)  IS  SET 
A  DEC'S  CONVENTION 
WAIT  A, ?B 

A  IS  THE  DEVICE  ADDRESS,  B  IS  AN 
ASSEMBLER  GENERATED  LABEL 

e#A 

B 


MAIN  PROGRAM 


TITLE  PASS1 

SBTTL  PASS1  OF  DATA  PROCESSING 
MAIN  PROGRAM  CONTROL  SECTION 
CSECT  MAIN 

ENTRY  POINT  NAME  IS  PASS1 
GLOBL  PASS1 
CALL  FOR  SYSTEM  MACROS 


.MCALL 

..V2..,. REGDEF 

.MCALL 

.PRINT, .FETCH,. EXIT 

.MCALL 

.ENTER,. CLOSE,. WRITE 

.MCALL 

.WRITW , .WAIT 

* 

. .V2. . 

* 

. REGDEF 

t 

• 

;  START 

PASS1 

MOV 

#ISR,g#KWIV 

• 

t 

THESE  ARE  ALL  THE  SYSTEM 
MACRO  CALLS 


THIS  IS  THE  FIRST  INTERRUPT  SERVICE 
ROUTINE  TO  TAKE  CARE  OF  PREFLIGHT 
DATA 

NO  FURTHER  INTERRUPT  IS  ALLOWED 
WHEN  THE  INTERRUPT  SERVICE  ROUTINE 


MOV 


#NOINT, 6IKWIV+2 
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Ifetch  #hndr,#device 

BCC  1$ 

MOV  #FERR,RO 

JMP  FAIL 

1$:  .ENTER  #AREA,#0,#FILE,#FILSIZ 


BCC  2$ 

MOV  #EERR,RO 

JMP  FAIL 

2$:  j  ERASE  THE  SCOPE  AND  GET 

MOV  #ERASE,@#VCST 

MOV  #VCRDY,§#VCST 

WAIT  VCST 

j  INITILIZE  SOME  PARAMETERS 

4lr  X 

CLR  RECNO 

CLR  WBUF 


;  IS  EXECUTING  (BR7) 

;  GET  THE  DK  HANDLER 
;  INTO  MEMORY 

;  SOMETHING'S  WRONG 
;  FATIAL  ERROR 

;  CREATE  PASS1.TMP  ON  THE  DISK 
;  ON  I/O  CHANNEL  0 
;  WITH  DEFAULT  FILE  SIZE 

;  CANNOT  CREATE  NEW  FILE 

:  FATIAL  ERROR 

READY 

;  ERASE  THE  603  SCREEN 
;  GET  READY,  Y  MODE,  FAST  INTENSIFY 
;  WAIT  TILL  IT  IS  READY 


X-COOR  OF  THE  LPSVC 
BLOCK  NUMBER  FOR  .WRITE  MACRO 
FLAG  TO  IDENTIFY  DOUBLE  BUFFER 
0  MEANS  BUFFER  1 
1  MEANS  BUFFER  2 


13$: 


tiov 

#KWENBL,0#KWST 

;  START  THE  CLOCK 

tiov 

#CHANO,@#ADST 

;  SAMPLE  ADC  CHANNEL  0 

WAIT 

ADST 

;  WAIT  FOR  ADC 

CMP 

g#ADBF, #ZEROV 

THE  ABOVE  NEEDS  SOME 

;  SEE  IF  IT  IS  A  VOLTAGE 

• 

• 

EXPLANATION 

• 

• 

THE  PROGRAM  WILL  EXECUTE  IN  IDLE 

• 

* 

• 

V 

UNTIL  THE  PUSH  BUTTON  AT  CHANNEL  0 

FIRES  A  PULSE 

• 

• 

THAT  MEANS  THE  FIRST 

WAVEFORM  APPEARS 

6le 

AND  DATA  ACQUISITION 
START 

PHASE  WILL  START 

OTHERWISE 

ENTER  AN  IDLE  LOOP  TO  KILL  SOME  TIME 
THE  IDLE  TIME  DEPENDS  ON  THE  PULSE  WIDTH 
OF  THE  PULSE  GENERATOR  AT  CHANNEL  0 
JUST  WANT  TO  MAKE  SURE  THAT  THE  SAME 
PULSE  IS  NOT  READ  TWICE 


MOV 

#ICNT1 , IN  1 

DEC 

INI 

BNE 

3$ 

;  LOOP 

• 

• 

f 

« 

NOW  SEE  IF  THE  SCREEN 
ERASE  IT  IF  NECESSARY 

IS  FULL 

4mp 

#4095., X 

BGE 

#^ASE,g#VCST 

;  NOT  YET 

MOV 

;  ERASE  SCREEN 

MOV 

#VCRDY,@#VCST 

;  SET  IT  UP  AGAIN 

WAIT 

VCST 

CLR 

X 

;  START  FROM  0  AGAIN 

BR 

13$ 

;  LOOP  AGAIN 

THE  REAL  THING  STARTS  HERE 


START: 


1$: 


dLR  S#KWST 

MOV  #KWSERV,§#KWIV 

MOV  #BUFSIZ,R2 

MOV  #BUFF 1 ,  R 1 


STOP  THE  CLOCK  FIRST 

PUT  A  NEW  INTERRUPT  SERVICE 

ROUTINE  TO  THE  LPSWK'S  VECTOR 

R2  IS  USED  AS  A  COUNTER  FOR  BUFFER 

R1  HAS  THE  ADDRESS  FOR  BUFFER  1 

START  FILLING  BUFFER  1  FIRST 

GET  BCD  TIME  FOR  THE  FIRST  WAVEFORM 


CALL  POLL 

;  IDLE  LOOP  TO  KILL  TIME 
MOV  0ICNT1 , IN  1 

DEC  INI 

BNE  1$ 

HERE  IS  THE  MAIN  LOOP 

THE  ROUTINE  WILL  LOOP  FOREVER  UNTIL  A  VOLTAGE  FALL  IS  DETECTED 
ON  CHANNEL  2,  THEN  IT  WILL  CLOSE  ALL  FILES  AND  EXIT 


LOOP: 

;  THIS 

IS  THE  MAIN  LOOP 

Aov 

#CHANO,§#ADST 

« 

SAMPLE  CHANNEL  0 

WAIT 

ADST 

CMP 

@#A DBF.IZEROV 

IS  IT  A  VOLTAGE  DROP 

BGT 

7$ 

• 

NO 

CALL 

POLL 

• 

OTHERWISE,  READ  IN  BCD  TIME 

7$: 

• 

• 

• 

THIS  IS  AN  IDEL  LOOP  TO 

KILL  TIME 

riov 

#ICNT2, IN  1 

4$: 

DEC 

INI 

BNE 

4$ 

MOV 

#CHAN2,@#ADST 

• 

GET  A  READING  FROM  CHANNEL  2 

WAIT 

ADST 

CMP 

@#ADBF , #ZEROV 

• 

IS  IT  A  DROP? 

BGT 

JMP 

ISt 

* 

• 

NO,  GO  ON 

YES,  THAT'S  ALL  FOLKS 

3$: 

CMP 

X,#4095. 

• 

IS  THE  SCREEN  FULL? 

BLE 

LOOP 

• 

NO 

MOV 

0ERASE,§#VCST 

• 

YES,  ERASE  THE  SCREEN 

MOV 

0VCRDY,@#VCST 

• 

SET  IT  UP  AGAIN 

WAIT 

VCST 

BR 

LOOP 

f 

GO  BACK  TO  LOOP 

♦ 

fH£  mAIn  Riufiiii  en6£  h£A£” 

SUBROUTINE  TO  READ  IN  TIME  CODE 

NOTE  THAT  THE  LPS  USES  NEGATIVE  LOGIC 

THEREFORE  ALL  BIT  PATTERNS  ARE  COMPLEMENTED 


SBTTL  POLL  TIME 
CSECT  POLL 
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POLL: 


ISR: 


KWSERV: 


TIMING  INFORMATION  COME  IN 
THRU  THE  DIO  PORTS 


4  BYTES  OF  BCD 


£lr 

MOV 

CALL 

MOV 

CLR 

MOV 


MOV 

BIS 

MOV 

CALL 

MOV 

BIS 

MOV 

CALL 

MOV 

BIS 

MOV 

CALL 

MOV 

BIS 

MOV 

CALL 


g#KWST 

#EOWF,R4 

STORE 

#ERASE,§#VCST 

X 

#VCRDY,g#VCST 


STOP  THE  CLOCK  FIRST 
PUT  OUT  A  NULL  FIRST 
STORE  THE  END  OF  WAVEFORM  MARK 
ERASE  THE  SCREEN  AS  WELL 
RESET  THE  X-COOR  COUNTER 
SET  THE  SCOPE  UP  AGAIN 


THE  DIGITAL  I/O  USES  NEGATIVE  LOGIC 

THEREFORE  PROGRAM  HAS  TO  SENT  OUT  THE  COMPLEMENTED  BIT  PATTERN 


@#DRIN,R4 

STORE 

#177775. @#DROUT 
#2,g#DRST 
@#DRIN , R4 
STORE 

#1 77773. @#DROUT 
#2,g#DRST 
@#DRIN,R4 
STORE 

#1 77767. §#DROUT 
#2,e#DRST 
@#DRIN,R4 
STORE 


A  -1  TO  PULL  IN  SECOND 
INPUT  THRU  DIGITAL  INPUT  PORT 
DATA  IN  R4 

STORE  IT  INTO  THE  BUFFER 
A  -2  TO  GET  MINUTES 


STORE  DATA  INTO  THE  BUFFER 
A  -4  TO  GET  HOUR 


;  A  -10  TO  GET  DAYS  OF  THE  YEAR 


MOV 

RETURN 


#KW  ENBL , 8  #KWST 


;  RESTART  THE  CLOCK 
;  GO  BACK 


THE  FIRST  INTERRUPT  SERVICE  ROUTINE 
TO  HANDLE  PRE -FLIGHT  DATA 

SBTTL  INTERRUPT  SERVICE  ROUTINES 
CSECT  ISR 

THIS  ROUTINE  JUST  DISPLAY  A  LINE  ON  THE  SCREEN 
WHENEVER  A  DATA  PULSE  COME  IN 
THIS  ROUTINE  TAKES  128  uS 
MOV  X.@#VCX 

MOV  #YVAL,g#VCY 

ADD  #XINC,X 

RTI  ;  JUST  TO  DISPLAY  A  POINT  ON  THE  SCOPE 

THIS  INTERRUPT  SERVICE  ROUTINE  HANDLES 
IN-FLIGHT  DATA 

IT  READS  THE  COUNT  FROM  LPSKW’S  BUFFER/PRESET 
REGISTER  AND  SAMPLE  CHANNEL  1  OF  THE  ADC 
THEN  DISPLAY  BOTH  OF  THEM  ON  THE  SCREEN  AND 
THEN  STORE  THEM  INTO  THE  CURRENT  BUFFER  IN  MEMORY 


MOV  @#KWBP,R4 

CALL  STORE 

MOV  ICHAN1 ,§#ADST 

WAIT  ADST 

MOV  g#ADBF,R4 

CALL  STORE 

MOV  X,g#VCX 


COUNT  IN  BUFFER/PRESET  REGISTER 

STORE  IT  INTO  THE  CURRENT  BUFFER 

SAMPLE  ADC  CHANNEL  1 

WAIT  TILL  IT  IS  READY 

GET  THE  READING 

STORE  IT  INTO  BUFFER 

LOAD  X-COOR 
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FAIL: 


EOT: 

2$: 

1$: 

6$: 

ill 

4$: 


MOV 

MOV 

ADD 

RTI 


#YVAL,0#VCY 

R4.##VCY 

ixInc.x 


S8TTL 

CSECT 


PRINT 

EXIT 


SBTTL 

CSECT 


fiLR 

TST 

BEQ 


MOV 

DEC 

BNE 

TST 

BEQ 

.WAIT 


LOAD  Y-COOR 

THE  ANALOG  CHANNEL  TOO 

INCREMENT  X 

RETURN  FROM  INTERRUPT 


ROUTINE  TO  HANDLE  FATIAL  ERRORS 

IT  WILL  PRINT  OUT  THE  ERROR  MESSAGE  AND  EXIT 

FAIL 

FAIL 

FATIAL  ERROR 


PRINT  ERROR  MESSAGE  AND  EXIT 


ERROR  MESSAGE  ADDRESS  IN  RO 
BYE.  BYE,  MISS  AMERICAN  PIE 


ENDING  ROUTINE 

IT  WILL  FILL  UP  THE  REST  OF  THE  CURRENT  BUFFER 
WITH  ASCII  "FI"  (FOR  FINISH) .  CLOSE  ALL  FILES 
AND  EXIT 

END  OF  TRANSMISSION 


EOT 

g#KWST 

R2 

1$ 

OTHERWISE 
FILL  WITH  "FI" 
EOJ, (R1 )+ 

R2 

wIuF 


IS 


;  WRITE 

.WRITW 

BCC 

MOV 

JMP 

BR 

.WAIT 

.WRITW 

BCS 

CLOSE 

PRINT 

EXIT 


SYSTEM  MACRO 
#AREA , #0 , IBUFF2 , 
5$  • 

IWERR.RO 

FAIL 

n 

IAREA , #0 , IBUFF 1 , 
6$ 

#0 

IENDMSG 


;  STOP  THE  CLOCK 

;  IS  THE  CURRENT  BUFFER  JUST  FULL? 

;  IF  YES,  IT  MUST  BE  A  BIG  COINCIDENT 

;  MOVE  END-OF-JOB  MARK 
;  WHICH  IS  "FI" 

;  UNTIL  THE  WHOLE  BUFFER  IS  FULL 
;  SEE  WHICH  BUFFER  THIS  IS? 

;  0  MEANS  IT  IS  BUFFER  1 
;  IN  CASE  THE  LAST  WRITE  IS  NOT  FINISH 
;  WAIT  FOR  IT  TO  COMPLETE 

IBUFSIZ, RECNO 

;  WRITE  OK? 

;  NO 

;  FATIAL  ERROR 
;  EXIT 

•  THIS  WRITES  OUT  BUFFER  1 
#BUFSIZ , RECNO 

;  ERROR 
;  CLOSE  FILE 
;  PRINT  ENDING  MESSAGE 
;  EXIT 


SBTTL  STORE  DATA 

STORE  DATA  POINTS  USING  DOUBLE  BUFFER  SCHEME 

DATA  TO  BE  STORED  IS  PASSED  IN  FROM  R4 

R1  IS  A  POINTER  TO  THE  CURRENT  AVALIABLE  LOCATION 

R2  CONTAINS  THE  NUMBER  OF  FREE  LOCATIONS  LEFT 

WBUF  IS  A  FLAG  TO  INDICATE  WHICH  BUFFER  IS  CURRENTLY 

BEING  USED  0  =  BUFFER  1.1*  BUFFER  2 

IF  NO  ROOM  IS  AVALIABLE  tN  THE  CURRENT  BUFFER, 

IT  WILL  SWITCH  BUFFER  AND  GUARANTEE  THE  DATA 
POINT  MUST  BE  SAVED 
.CSECT  STORE 
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STORE: 


TST 

BNE 


R2 

3$ 


ANY  ROOM  LEFT 

YES,  EVERYTHING  IS  OK 


OTHERWIES  !!!M 


1$: 


4$: 

2$: 

3$: 


£lr 


TST 

BNE 

WAIT 

WRITE 


BCS 

MOV 

INC 

BR 

.WAIT 

.WRITE 

BCS 

MOV 

CLR 

BR 

MOV 

JMP 

ADD 

MOV 

MOV 

MOV 

DEC 

RETURN 


@#KWST 


ONE  BUFFER  IS  FULL 

HAS  TO  WRITE  IT  OUT  AND  SWITCH  BUFFER 
*•*  STOP  THE  CLOCK  ••• 

SEE  WHICH  BUFFER  IS  FULL 
1  MEANS  BUFFER  2 

IN  CASE  THE  LAST  WRITE  IS  NOT  FINISH 


WBUF 
1$ 

#0  .  x., 

#AREA , #0 . IBUFF 1 , #BUFSIZ , ftECNO 
ISSUE  AN  ASYNCHRONOUS  WRITE 
CONTROL  WILL  PASS  BACK  TO  PROGRAM  RIGHT  AFTER 
THE  WRITE  REQUEST  IS  QUEUED 


,1 


ERROR 

RESET  THE  POINTER  TO  BUFFER  2 

SET  FLAG  TO  INDICATE  BUFFER  2  IS  IN  USE 


UFF2.R1 
WBUF 
2 $ 

#0  :  THIS  TIME  IS  TO  WRITE  OUT  BUFFER  2 

#AREA , #0 , IBUFF2 , #BUFSI Z , RECNO 
4$  ;  ERROR 

#BUFF 1 , R 1 
WBUF 
2$ 

#WERR , RO 
FAIL 

0RECLEN, RECNO 
#BUFSIZ,R2 
IKWENBL , 6#KWST 
R4, (R 1 )+ 

R2 


RESET  POINTER  TO  BUFFER  1 
INDICATE  BUFFER  1  IS  IN  USE 


WRITE  ERROR,  FATIAL 
UPDATE  BLOCK  COUNTER 
RESET  POINTS  COUNTER 
RESTART  THE  CLOCK 
STORE  POINT  INTO  BUFFER 


DEVICE: 

FILE: 

AREA: 

RECNO: 

WBUF: 

EOJ: 

INI: 

X: 


FERR: 

EERR: 

WERR: 

ENDMSG: 

BUFFI: 

BUFF2: 

HNDR: 


.SBTTL  DATA  DEFINITIONS 
. CSECT  DATA 
• 

,’nlist 

. RAD50/DK  / 

. RAD50/DK  PASS1  TMP/ 
.BLKW  10. 

.WORD  0 
.BLKW  1 
.WORD  "FI 
.WORD  0 
.WORD  0 

;  ERROR  MESSAGES 

!asciz/  NO  DEVICE  / 

. ASCIZ/  CANNOT  CREATE  / 

. ASCIZ/  WRITE  ERROR  / 
.ASCIZ/  END  OF  PASS1  / 


.EVEN 
.BLKW 
.BLKW 
.*2 
.LIST 

.END  PASS1 


4096. 

4096. 


DEVICE  NAME 

FILE  NAME 

AREA  FOR  WRITE 

WRITE  RECORD  NUMBER 

FALG  TO  INDICATE  WHICH  BUFFER 

END  OF  JOB  MARK 

FREE  SLOT  FOR  COUNTING 

X-COOR 


FETCH  ERROR 
ENTER  ERROR 
WRITE  ERROR 
ENDING  MESSAGE 

BUFFER  1,  4K 
BUFFER  2,  4K 
PLACE  TO  PUT  DK  HANDLER 


oooooooooo 
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102 


.ENDM 
.MACRO 
RTS 
.ENDM 
•  LIST 


RETURN 

17 


SOME  CONSTANT  DEFINITIONS 


;  RETURN  FROM  SUBROUTINE 


£rRWD  =  52 
INSIZ  =  4096. 
OUTSIZ  =  4096. 
F180  =  560. 
FEQLIM  =  ‘ 
SLPLIM  = 

;  LPS  ADDRESSES 
ADBF  =  170402 
INCNT  =  16. 


HEADLN  =  256. 


SYSTEM  ERROR  WORD 

INPUT  BUFFER  SIZE  =  4K 

OUTPUT  BUFFER  SIZE  =  4K 

COUNT  CORROSPONDS  TO  180  HZ 

MORE  THAN  8  PTS  >  180  HZ 

MORE  THAN  8  PTS  HAVE  NON-RISING  SLOPE 

ADC  DATA  REGISTER 
INPUT  BUFFER  BLOCK  LENGTH 
4K  =  16.  BLOCKS 
HEADER  =  128  WORD  LONG 
=  256.  BYTES 


CANNOT  USE  THE  NAME  "PASS2"  BECAUSE  THE  CALLING  FORTRAN 
PROGRAM  REQUIRES  THIS  NAME 

OTHERWISE  WILL  GET  MULTIPLE  DEFINITION  ERROR  FROM  LINK 
TITLE  PASS20 

SBTTL  PASS 2  OF  DATA  PROCESSING 

THIS  PASS  READS  IN  DATA  OBTAINED  FROM  PASS1 
INTO  THE  INPUT  BUFFER.  DECODES  4  WORDSIONLY  THE  LOWER  BYTE 
OF  EACH  WORD)  OF  TIMING  INFORMATION  INTO  BINARY  FORM 
SPLIT  THE  INCOMING  DATA  INTO  INDIVIDUAL  WAVEFORMS 
AND  THROWS  OUT  UNNECESSARY  DATA  (  DATA  THAT  HAS  A  MAJORITY 
OF  NON-RISING  SLOPES  AND  A  MAJORITY  OF  FREQUENCIES  >  180  HZ) 
FINALLY  IT  ROUNDS  THE  DATA  INTO  THE  NEXT  256  WORD  BOUNDARY 
AND  WRITE  THE  WAVEFORM  TO  THE  OUTPUT  FILE  "PASS2.TMP" 


SYSTEM  MACRO  CALLS 


MCALL  . . V2. . , . REGDEF 
MCALL  .FETCH,. ENTER,. LOOKUP 
MCALL  . READW, .WRITW 
MCALL  .EXIT,. CLOSE,. PRINT 
ENRTY  POINT 
GLOBL  PASS 2 


.  V2. . 
REGDEF 


CONTROL  SECTION  DEFINITION 
CSECT  PASS20 

GET  RK05  HANDLER  IN  MEMORY 
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PASS2:  .FETCH  IHNDR , IDEVICE 

BCC  1| 

MOV  #FERR,RO 

JMP  FAIL 

1$:  .LOOKUP  IINAREA , #0 , IINFIL 

BCC  2% 

TSTB  £#£RRWD 

BEQ  3$ 

MOV  INERR.RO 

JMP  FAIL 

3$:  MOV  #AERR,RO 

JMP  FAIL 

2$:  .ENTER  IOAREA ,  #  1 , IOUTF IL , #- 1 


IEERR.RO 

FAIL 


FETCH  OK? 

NO.  FATIAL  ERROR 
FETCH  FAIL 

OPEN  'PASS1 .TMP'  ON  CHANNEL  0 
FOR  INPUT 
IS  IT  OK? 

WHAT'S  WRONG? 

0  MEANS  CHANNEL  ACTIVE 
NO  FILE 

CHANNEL  IS  ACTIVE 

CREATE  OUTPUT  PTLE 
"PASS2.TMP"  ON  DISK 
ENTER  OK? 

CREATE  ERROR 


INITILIZE  SOME  PARAMETERS 


.SBTTL  INIT  PARAM 


INIT  PARAMETERS 


LR  INBLK  ,*  INPUT  BLOCK  COUNTER 

LR  OUTBLK  ;  OUTPUT  BLOCK  COUNTER 

OV  #INBUF, R 1  ;  R1  HAS  THE  ADDRESS  OF  THE 

;  INPUT  BUFFER 

OV  #INSIZ,R2  ;  R 2  HAS  THE  SIZE  OF  THE  INPUT  BUFFER 

BUT  THE  OUTPUT  BUFFER  IS  ACTUALLY  4K  -  128  WORDS  (HEADER) 

PUT  THIS  NUMBER  ON  "OUTCNT" 

OV  #OUTSIZ  -  128.. OUTCNT  ;  UK  -  128.  WORDS 

LR  ENDFLG  ;  END  FILE  FLAG 

SBTTL  PROCESS  HEADER 

THIS  IS  DONE  BY  THE  FORTRAN  SUBROUTINE  "QUEST" 


#OUTBUF.R3 


R3  POINTS  TO  BEGINNING  OF  OUTBUF 


nwi  rvwiww t  •  it  j  ■  j  i  isahav/  *  v/  vi  vu  x  uv 

;  BEGIN  AT  WHERE  A  ZERO  IS  INSERTED 

;  SKIP  THE  HEADER  UNTIL  THE  POSITION  IN  WHICH  THE  "0" 

;  IS  PLACED 

Add  #236. ,R3 

CLR  (R3)+  5  PUT  A  ZERO  IN  THERE 

MOV  R3.BADDR  ;SAVE  THIS  ADDRESS  FOR  LATER  USE 

J SBTTL  READ  IN  THE  1'ST  BUFFER 

;  READ  IN  THE  FIRST  BUFFER  JUST  TO  GET  STARTED 

.READW  #INAREA,#0,R1,R2, INBLK 

BCC  6$  ;  READ  ERROR  ? 

MOV  IRERR.RO 


irerr.ro 

FAIL 

IINCNT, INBLK 


;  FATIAL 
;  UPDATE  INBLK 


THE  REAL  PROCESSING  STARTS  HERE 


SBTTL  START 


;  START  PROCESSING 


J 

i . 


104 


LOOP: 


.CSECT  START 


NOREAD:  TST 


lSop 

R4,(R3)+ 

ZFLAG 

READ 

DECODE 

R4.SEC 

READ 

DECODE 

R4 


SAVE  IT  IN  "SEC" 


R4 
4$ 

ADD  #3600., SEC  ;  YES,  ADD  3600  SECONDS 

MOV  SEC.R4  ;  NOW  "SEC"  HAS  THE  BINARY  SECONDS 

;  OF  THE  BCD  TIME  CODE 

MOV  R4. (R3)  ;  PUT  IT  INTO  THE  OUTPUT  BUFFER 

ADD  #16., S3  ,*  SKIP  TO  THE  BEGINNING  POSITION 

:  FOR  HOLDING  DATA  POINTS 

;  LED  ROUTINE  WILL  DISPLAY  THE  CONTENT  OF  R4  ON  THE  LED  DISPLAY 
:  ON  THE  LPS  IN  DECIMAL  FORM 
CALL  LED 

:  READ  IN  THE  DAYS  OF  YEAR,  THIS  IS  NOT  USED  AT  PRESENT 
CALL  READ 


READ  IN  THE  1'ST  POINT 

RETURNS  IN  R4 

UNTIL  IT  FINDS  A  ZERO 

WHICH  IS  THE  END-OF -WAVEFORM  MARK 

OTHERWISE  JUST  KEEP  LOOPING 
PUT  IT  IN  THE  OUTPUT  BUFFER  FIRST 
CLEAR  THIS  FLAG 

THE  t 'ST  TIMING  WORD  RETURNS  IN  R4 
ONLY  THE  LOWER  BYTE  IS  USED 
DECODE  IT  INTO  BINARY  FORM 
SAVE  IT 
NEXT 

IS  IT  ZERO? 

CONVERT  MINUTES  TO  SECONDS 


NEXT 

SINCE  LPS  IS  IN  NEGATIVE  LOGIC 
THEREFORE  COMPLENT  IT 
IS  IT  SET? 


THE  FOLLOWING  SECTION  OF  CODE  IS  FOR  REFINING  THE  DATA  POINTS 
IT  READS  IN  16.  PAIRS  OF  POINTS  AT  A  TIME 
IF  THIS  16.  PAIRS  OF  POINTS  SATISFY: 

A)  HAVE  A  MAJORITY  OF  POINTS  >  180  HZ  AND 

B)  HAVE  A  MAJORITY  OF  NON-RISING  SLOPES 

THEN  THE  REST  OF  THE  DATA  POINTS  IN  THE  WAVEFORM  IS 
VERY  LIKELY  TO  BE  USELESS  FOR  CALCULATION  BECAUSE 
CONDUCTIVITY  CALCULATIONS  INVOLVES  ONLY  THE  RISING 
SLOPES  OF  THE  WAVEFORMS 

THEREFORE  THE  REST  OF  THE  DATA  POINTS  WILL  BE  DISCARDED 
TO  SAVE  STORAGE  SPACE 


READ10:  dLR  FEQFLG 


SLPFLG 

#31., TSTCNT 

READ 

R4 


THIS  FLAG  IS  USED  TO  INDICATE  HOW  MANY 

POINTS  ARE  ABOVE  180  HZ 

THIS  FLAG  IS  USED  TO  INDICATE  HOW  MANY 

SLOPES  ARE  NON-RISING 

31  POINTS  LEFT 

30  LEFT  (15  PAIRS) 


IF  IT  EQUALS  TO  0  MEANS  THATS  THE  END  OF  THE  WAVEFORM 
NOTHING  HAS  BEEN  THROWN  OUT 


f 
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SAVE  10: 


BEQ  SAVE  10  ;  TAKE  CARE  OF  THE  ROUNDING 

MOV  R4,(R3)+  ;  SAVE  IT  IN  THE  OUTPUT  BUFFER 

DEC  OUTCNT  ;  UPDATE  THE  COUNTER 

;  HERE  IS  A  LITTLE  TRICK  TO  DETERMINE  IF  THE  SLOPE  IS  RISING 
;  OR  FALLING 

;  SINCE  THE  COUNT  IS  INVERSELY  PROPORTIONAL 
;  TAKE  THE  DIFFERENCE  BETWEEN  2  CONSECUTIVE  COUNT 
;  DEPENDING  ON  WETHER  THE  DIFFERENCE  IS  POSITIVE  OR  NEGATIVE 
;  CAN  DETERMINE  THE  SLOPE  IS  RISING  OR  FALLING 


R4.R5 

#F180,R4 


;  HAS  TO  GET  THE  1'ST  POINT  IN 
;  SEE  IF  IT  IS  >  180  HZ 

:  IF  >  INCREMENT  FUG 


v  a  |  s  inviibiibni  r  i 

NOW  ENTER  THE  LOOP  FOR  THE  REST  OF  THE  31  POINTS 


U.L  READ 

3V  R4,(R3)+ 
iC  OUTCNT 

iC  TSTCNT 

ILL  READ 

3T  R4 

:q  save io 

3V  R4,(R3) 

SC  OUTCNT 

COMPARE  FREQUENCY  AND  SLOPE 


READ  IN  R4 
SAVE  IN  OUTBUF 

UPDATE  BOTH  COUNTERS 


ANALOG  TOO 


CMP  #F180,R4 

BLT  6$ 

INC  FEQFLG 

SUB  R5.R4  ;  IF  LAST  COUNT  >  PRESENT  COUNT 

;  =>  RISING  SLOPE 

BLT  7$ 

INC  SLPFLG 

MOV  (R3)+,R5 

DEC  TSTCNT 

BNE  5$ 

;  SEE  IF  BOTH  LIMITS  EXCEED 

:  IF  THEY  DO.  THE  REST  OF  THE  POINTS  CAN  BE  THROWN  OUT 

CMP  #FEQLIM, FEQFLG 

BGE  READ 10 

CMP  #SLPLIM, SLPFLG 

BGE  READ 10 

;  IF  EXECUTION  REACHES  HERE.  BOTH  LIMITS  HAVE  BEEN  EXCEEDED 
;  THE  ZERO  FLAG  IS  USED  TO  INDICATE  WHETHER  IT  NEEDS  A  READ  OR  NOT 
;  IF  THE  n0"  HAS  ALREADY  BEEN  READ,  THERE  IS  NO  NEED  TO  READ  AGAIN 
;  THIS  IS  TO  KEEP  THE  PROPER  ALIGNMENT  OF  DATA 
INC  ZFLAG  :  BACK  SPACE  0  OR  NOT 

;  THE  FOLLOWING  SECTION  OF  CODE  WRAPS  UP  THE  PIECES 
;  IT  ROUNDS  UP  THE  REST  OF  THE  BUFFER  TO  THE  NEXT 
;  256  WORDS  BOUNDARY  AND  WRITE  IT  TO  THE  DISK 

MOV  #OUTSIZ,R5  ;  DETERMINE  HOW  MANY  POINTS  HAVE 

;  BEEN  ENTERED 

SUB  OUTCNT, R5 

MOV  R5.R4 

SUB  #128. ,R4 

:  DIVIDED  BY  2  GIVES  NUMBER  OF  POINT  PAIRS 

ASR  R4  ;  THIS  IS  THE  NUMBER  OF  PAIRS  OF  POINTS 

MOV  R4.NPTS  ;  SAVE  IT 

CLR  R4 

INC  R4  ;  CALCULATE  TO  THE  NEXT  256  WORDS 

SUB  #256., R5 

BGT  1$ 

TST  R5 

BEQ  2$  *,  OR  MAY  BE  IT  IS  JUST  RIGHT  AT 
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1  i ; 

■  t  ; 


256  WORDS  BOUNDARY  ?? 
CLEAR  THE  REST 


CLR  (R3)+  5  CLEAR  THE  REST 

INC  R5 

BNE  3$ 

MOV  R 4, NBLK  ;  SAVE  THE  NUMBER  OF  BLOCKS 

CLR  R5 

;  THE  SYSTEM  MACRO  WRITE  REQUIRES  THE  BUFFER  SIZE  IN  WORDS 
:  THEREFORE  HAS  TO  CALCULATE  256  •  NUMBER  OF  BLOCKS 
ADD  #256., R5  ;  GET  BACK  NBLKS»256  TO  WRITE 

DEC  R4 

BNE  4$ 

;  WRITE  (NUMBER  OF  BLOCKS)  AND  (NUMBER  OF  DATA  POINT  PAIRS) 

;  TO  THE  HEADER 
MOV  #OUTBUF.R4 

ADD  #248., R 4 

MOV  NPTS, (R4) 

ADD  #4,R4 

MOV  NBLK, (R4) 

;  WRITE  OUT  THE  WAVEFORM  TO  DISK 
.WRITW  #OAREA , # 1 , #OUTBUF , R5 , OUTBLK 


VVAHbA  a  If  I  a  A  UUl  •  11. 

NBLK, OUTBLK 

#OUTSIZ-128.,OUTCNT 

BADDR.R3 

ZFLAG 

lSop 

R4 

NOREAD 


UPDATE  COUNTER 
RESET  OUTPUT  SIZE  COUNT 
RESET  BEGINNING  ADDRESS 
SEE  IF  I  NEED  TO  BACK  SPACE  0 


'm’msymNmmt' 


THIS  ROUTINE  WILL  RETURN  A  DATA  POINT  IN  R4 
IT  AUTOMATICALLY  MAINTAINS  INPUT  BUFFERING 
A  4K  INPUT  BUFFER  IS  USED 
IF  DATA  FROM  INPUT  BUFFER  IS  EXAUSTED 
IT  WILL  READ  IN  ANOTHER  BUFFER  FROM  THE  FILE  "PASS1.TMP" 
UNTIL  END  OF  FILE  IS  REACHED 
SBTTL  READ 


IS  THIS  THE  LAST  BUFFER  OF  THE  FILE 


.CSECT  READ 

TST  R2  ;  R2  HAS  INPUT  BUFFER  COUNT 

;  IS  THERE  ANYTHING  LEFT  ? 
fiNE  1$  ;  YES 

TST  ENDFLG  ;  IS  THIS  THE  LAST  Bl 

BEQ  2$ 

JMP  EOF  :  YES,  LAST  BUFFER 

;  OTHERWISE.  READ  IN  ANOTHER  BUFFER  FROM  "PASS1.TMP" 

MOV  #INBUF,R1 

. READW  #INAREA,#0,R1,#INSIZ,INBLK 
BCC  3$  ;  READ  OK? 

TSTB  elERRWD  ;  WHAT'S  WRONG? 

MO?  I^ERR.RO  ;  READ  ERROR 

JMP  FAIL  ;  READ  ERROR,  FATIAL 

TST  RO 

BNE  5$ 

JMP  EOF 
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R0,R2 

ENDFLG 

#INCNT, INBLK 

#INSIZ , R2 

(R1)+,R4 

R2 

R4.EOJ 

6$ 

EOF 


END  OF  FILE,  RO  CONTAINS  ACTUAL 

NO  OF  POINTS  READ 

SET  END-OF-FILE  FLAG 

GO  ON  WITH  THE  LAST  BUFFER 

UPDATE  INPUT  BLOCK  COUNTER 

4K  =  16  BLOCKS 

RESET  COUNTER 

RETURN  DATA  IN  R4 

UPDATE  COUNTER 

IS  IT  THE  END  OF  FLIGHT  MARK 
YES,  THE  END 


RETURN 


DECODE  BCD  TIME 
SBTTL  DECODE 

DECODE  THE  LOWER  BYTE  OF  R4  WHICH  CONTAINS  2 
BCD  DIGITS  AND  RETURNS  THE  BINARY  RESULT  IS  R4 
THE  BCD  NUMBER  IS  IN  NEGATIVE  LOGIC 
THEREFORE  HAS  TO  BE  COMPLEMENTED  FIRST 
CSECT  DECODE 

THIS  MIGHT  LOOK  CLUMSY  BUT  IT  IS  THE  MOST  EFFICENT  WAY 
BELIEVE  IT  OR  NOT 


DECODE: 

£lr 

R5 

COM 

R4 

RORB 

R4 

BCC 

1$ 

INC 

R5 

1$: 

RORB 

R4 

BCC 

2$ 

ADD 

#2,R5 

2$: 

RORB 

BCC 

R4 

ADD 

I4.R5 

3$: 

RORB 

R4 

BCC 

4$ 

ADD 

#8.,R5 

4$: 

RORB 

BCC 

R4 

5$ 

t  5$: 

ADD 

#10. ,R5 

RORB 

R4 

BCC 

6$ 

ADD 

#20., R5 

6$: 

i 

RORB 

R4 

BCC 

7$ 

1 

ADD 

#40., R5 

7$: 

RORB 

R4 

BCC 

8$ 

i 

ADD 

#80.,R5 

j  8$: 

MOV 

R5,R4 

RETURN 

R5  IS  TO  HOLD  THE  RESULT  TEMPORARY 

R4  HAS  THE  BCD  NUMBER 

R5  WILL  HAVE  THE  BINARY  NUMBER 


; RETURN  RESULT  IN  R4 


;  CLOSE  ALL  FILES,  PRINT  OUT  ENDING  MESSAGE  AND  EXIT 
.SBTTL  EOF  ;  END-OF-FILE 


EOF: 


.CLOSE  #0 


;  CLOSE  INPUT  FILE 


108 


.CLOSE  #1 
.PRINT  IENDMSG 
.EXIT 


CLOSE  OUTPUT  FILE 
END  MESSAGE 


FAIL: 


LED: 


SI 


. SBTTL  FATIAL  ERROR 
;  PRINT  OUT  ERROR  MESSAGE  AND  EXIT 

.’print 

•  EXIT 


MESSAGE  ADDRESS  IN  RO 


.SBTTL  LED 
.CSECT  LED 

;  R4  CONTAINS  THE  BINARY  NUMBER  TO  BE  DISPLAYED 

;  DIVIDE  IT  INTO  POSITIONAL  FORM (DECIMAL)  AND  PUT  IT  OUT  ON  THE  LED 
:  DISPLAY  ON  THE  LPS 
CLR  LED1 

;  CLEAR  ALL  6  DIGITS  FIRST 
MOV  #17,@#ADBF 

#417,@#ADBF 
#1017,@#ADBF 
#14l7,f#ADBF 
#2017,e#ADBF 
#2417,8#ADBF 
#-1 .LED2 
LED2 

#10. , R4  ;  GET  REMAINDER 

2$ 

#10.. R4 

UP  A  WHOLE  WORD  AND  GO 
R4.LED3 
LED1.LED3+1 
-  DBF 


MOV 
MOV 
MOV 
MOV 
MOV 
MOV 
INC 
SUB 
BGE 
ADD 
;  SET 
MOVB 
MOVB 
MOV 
INC 
MOV 
TST 
BNE 
RETURN 
:  END  LED 


LED3,e*Al 

LED1 

LED2.RH 

R4 

1$ 


ILLUMINATE 


LED1: 

.WORD 

0 

LED2: 

.WORD 

0 

LED3: 

.WORD 

0 

! SBTTL 

DATA 

.CSECT 

DATA 

EOJ: 

.WORD 

"FI 

INBLK: 

.WORD 

0 

OUTBLK : 

.WORD 

0 

OUTCNT : 

.WORD 

0 

TSTCNT : 

.WORD 

0 

BADDR : 

.WORD 

0 

SEC: 

.WORD 

0 

NBLK: 

•  WORD 

0 

NPTS: 

•  WORD 

0 

ZFLAG : 

.WORD 

0 

ENDFLG : 

•  WORD 

0 

SLPFLG : 

.WORD 

0 

FEQFLG: 

.WORD 

0 

INAREA : 

.BLKW 

10. 

OAREA: 

.BLKW 

10. 

INBUF : 

.BLKW 

4096 

FINISH  MARK 
INPUT  BLOCK  COUNT 
OUTPUT  BLOCK  COUNT 
OUTPUT  BUFFER  WORD  COUNT 
TEST  COUNT  FOR  31  POINTS 
TO  HOLD  BEGIN  ADr 
STORE  SECOND 
NUMBER  OF  BLOCKS 
NUMBER  OF  POINTS 
FREQ/SLOPE  FLAG 
END-OF-FILE  FUG 
SLOPE  FUG 
FREQUENCY  FUG 
SYSTEM  MACRO  AREA 
FOR  OUTPUT 
INPUT  BUFFER  4K 


DRESS  FOR  OUTBUF 


NAME  THE  OUTPUT  BUFFER  BY  CSECT  NAME  OBUF 
TO  MAKE  IT  KNOWN  TO  FORTRAN  SUBROUTINE  QUESTN 
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.CSECT  OBUF 
OUTBUF:  .BLKW  4096. 

. NLIST 
.CSECT  DATA 

INFIL:  .RAD50/DK  PASS1  TMP/ 
OUTFIL:  .RAD50/DK  PASS2  TMP/ 
DEVICE:  . RAD50/DK  / 

FERR:  .ASCIZ/NO  DEV/ 

NERR:  .ASCIZ/NO  FILE/ 

AERR:  .ASCIZ/CHAN  ACT/ 

EERR:  .ASCIZ/CANT  CREATE/ 

RERR:  .ASCIZ/READ  ERR/ 

OERR:  .ASCIZ/NO  ROOM/ 

WERR:  .ASCIZ/WRITE  ERR/ 

ENDMSG:  .ASCIZ/END  PASS2/ 
.EVEN 

HNDR:  .+2 

.LIST 


;  OUTPUT  BUFFER  4K 


INPUT  FILE  NAME 
OUTPUT  FILE  NAME 

FETCH  ERROR 
NO  FILE 
CHANNEL  ACTIVE 
ENTER  ERROR 
READ  ERROR 
OUTPUT  BUFFER  FULL 
WRITE  ERROR 
END  MESSAGE 


.END  PASS2 
SUBROUTINE  QUESTN 
DECURE  COMMON  AREA 

COMMON  /OBUF/DATE  (  2 ) ,  ST YPE ,  SNUMB ,  USITE ,  RF ,  RCAL ,  R 1 ,  R2 ,  L , DFDTCL , 
1DTSW , DVSW , VSWN . VSWP , DUM 1 , IDEN , TVO , ALT , VERVEL , ISAT , DUM2 , 

2SIG (S , 4 ) , DUM3 . DUM4 , DUM5 , ZERO , TIME (2 ) , NPTS , NBLKS 
REAL  lAsiTE.L.ISAT 
INTEGER*2  l6EN(8) 

INTEGER*2  STYPE, SNUMB 
INTEGER*4  NPTS. NBLKS 

SUBROUTINE  TO  FILL  INFORMATION  BLOCK  FOR  HEADER  OF  EACH  RECORD 


C 


100 

WRITE ( IOUT ,  1) 

1  FORMAT ( '  PROGRAM  TO  INPUT  HEADER  INFORMATION  FOR  EACH  FLIGHT*) 
C 

C  GET  DATE 
C 

WRITE (IOUT.  2) 

2  FORMAT ( *  DATE— DDMMMYY ' ) 

READ(IN , 3)  DATE(1),DATE(2) 

3  FORMAT (2A4) 

C 

C  GET  SENSOR  TYPE 
C 

WRITE (IOUT,  4) 

4  FORMAT ( *  SENSOR  TYPE  BP  FOR  BLUNT  PROBE —  GC  FOR  GERDIEN ' ) 
READ(IN,5)  STYPE 

5  FORMAT (A2) 

C 

C  GET  SENSOR  NUMBER 
C 

WRITEUOUT.  6) 

6  FORMAT ('  SliNSOR  NUMBER--  13') 

READ(IN.T)  SNUMB 

7  FORMAT (t 3) 

C  GET  LAUNCH  SITE 


IN  =  5 
IOUT  =  7 
ALT  =  -1.0 
DO  100  I  =  1,8 
DO  100  J  =  1,4 
SIG(I.J)  =  0.0 


ono  ooo  ooo  ooo  non  ooo  ooo  ooo  ooo 
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WRITE(I0UT,8) 

8  FORMAT ( '  LAUNCH  SITE  —  XXXX'} 

READ(IN,9)  LASITE 

9  FORMAT (A4) 

GET  FEEDBACK  RESISTOR 
WRITE (IOUT, 10) 

10  FORMAT ( ’  FEEDBACK  RESISTOR  SIZE  -RF-  F7.2') 

READ(IN.II)  RF 

11  FORMAT (E7. 2) 

GET  CALIBRATION  RESISTOR 
WRITE(IOUT, 12) 

12  FORMAT ('  CALIBRATION  RESISTOR  SIZE  -RCAL-F7.2’) 

READ(IN,13)  RCAL 

13  FORMAT (E7.2) 

COLLECTOR  RADIUS 
WRITECI0UT.14) 

14  FORMAT ( '  COLLECTOR  RADIUS-R  FOR  BP— RI  FOR  GC  F4.1') 

READCIN , 15)  RI 

15  F0RMATCF4. 1 ) 

GUARD  RADIUS 

WRITE(IOUTt 16) 

16  FORMAT ( '  GUARD  OR  OUTER  PLATE  RADIUS-  R  FOR  BP-RO  FOR  GC  F4.1') 
READ(IN,17)  R2 

17  FORMAT (F4. 1 ) 

ELECTRODE  LENGTH 
WRITEU0UT.18) 

18  FORMATC  ELECTRODE  LENGTH  IN  CM  -  ENTER  0  FOR  BP  F4.11) 

READ(IN . 19)  L 

19  FORMATCr 4. 1 ) 

DF/DT  CAL 

WRITE (IOUT, 20) 

20  FORMATC  DF/DT  CAL  F6.2') 

READ(IN,21)  DFDTCL 

21  F0RMAT(F6.2) 

DELTA  TIME  SWEEP 
WRITE (IOUT, 22) 

22  FORMATC  DELTA  TIME  SWEEP  F6.2’) 

READ(IN,23)  DTSW 

23  FORMAT(F6.2) 

DELTA  VOLTAGE  SWEEP 
WRITE (IOUT, 24) 

24  FORMATC  DELTA  VOLTAGE  SWEEP  F6.2') 

READ (IN. 25)  DVSW 

25  F0RMAT(F6.2) 

NEGATIVE  MAXIMUM  VALUE  OF  VOLTAGE  SWEEP 


WRITE (IOUT, 26) 


nononnnnonnnnonnnnnnonnonnnnoonnnnonnnnnnn  oo  oo  ooo 
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26  FORMAT ( '  NEGATIVE  MAXIMUM  VALUE  OF  SWEEP  VOLTAGE  F6.2') 

READ(IN,27)  VSWN 

27  FORMAT (F6. 2) 

POSITIVE  MAXIMUM  VALUE  OF  VOLTAGE  SWEEP 
WRITE (IOUT, 28) 

28  FORMAT ( '  POSITIVE  MAXIMUM  VALUE  OF  SWEEP  VOLTAGE  F6.2') 

READ(IN,29)  VSWP 

29  FORMAT (F6. 2) 

SPECIAL  IDENTIFICATION  OR  INFORMATION 
WRITE (IOUT, 30) 

30  FORMAT ( '  SPECIAL  IDENTIFICATION  OR  INFORMATION  8A2') 

READ(IN, 31 )  (IDEN(I9),I9s1,8) 

31  FORMAT (8A2) 

END  OF  QUESTION  AND  ANSWER  SUBROUTINE 
WRITE(I0UT,32) 

32  FORMAT ( '  END  OF  QUESTION  AND  ANSWER  SECTION  OF  PROGRAM') 

RETURN 

END 

MAIN  PROGRAM  FOR  PASS3  OF  DATA  PROCESSING 
PROGRAM  LAY  OUT: 

THIS  PROGRAM  WILL  CALL  THE  FOLLOWING  SUBROUTINES 

SUBROUTINE  INIT  -  INITILIZE  ALL  THE  PROGRAM  PARAMETERS  INCLUDING 
I/O  LOGICAL  UNIT  NUMBERS  AND  CONSTANTS,  OPEN  INPUT 
AND  OUTPUT  FILES. 

SUBROUTINE  LED  -  DISPLAY  TIMING  INFORMATION  ON  THE  LED  DISPLAY 
OF  THE  LPS  UNIT. 

SUBROUTINE  READIN  -  READ  IN  THE  NEXT  WAVEFORM.  SINCE  WAVEFORMS 

ARE  IN  VARIABLE  LENGTH,  THEREFORE  IT  HAS  TO  READ  IN  THE  1 'ST 
256  WORD,  FIND  OUT  THE  LENGTH  OF  THE  WAVEFORM  FROM  THE  HEADER 

SUBROUTINE  DISPR  -  DISPLAY  DATA  ON  THE  TEKTRONIX  603  SCOPE 

ACCORDING  TO  THE  FREQUENCY  LIMITS.  A  FLAG  IS  USED  TO 
INDICATE  WHEATHER  THE  ANALOG  CHANNEL  IS  TO  BE  DISPLAIED  OR  NOT. 

SUBROUTINE  EXPAND  -  DETERMINE  WHEATHER  THE  WAVEFORM  IS  TO  BE  EXPANDED 
FOR  DIFFERENT  LIMITS. 


SUBROUTINE  SLOPE  -  GET  THE  BREAK  POINTS  OF  THE  WAVEFORM  FROM  THE 

TERMINAL  AND  PERFORM  A  WEIGHTED  LEAST  SQUARE  STRAIGHT  LINE 
FIT  TO  THE  SLOPE  OF  THE  WAVEFORM,  THEN  DISPLAYS  THE  FITTED 
STRAIGHT  LINE  ON  THE  SCREEN 


SUBROUTINE  ENDW  -  END  OF  WAVEFORM  ROUTINE.  CALCULATE  THE  CONDUCTIVITY 
VALUE  FROM  THE  SLOPE,  RESET  THE  PARAMETERS.  WRITE  THE  HEADER 
BACK  TO  THE  INPUT  FILE  AND  WRITE  CONDUCTIVITY  INFO  TO  OUTPUT  FILE 


DATA  BASE  DEFINITIONS 


COMMON  DATA  -  INTEGER  ARRAY  IDATA(4096)  WHICH 

IS  THE  INPUT  BUFFER  TO  HOLD  ALL  THE  DATA  POINTS 
COMMON  PARM  -  CONTAIN  ALL  THE  PARAMETERS  AND  CONSTANTS 
I IN  -  LOGICAL  UNIT  NUMBER  FOR  INPUT  FILE  (19) 

IOUT  -  LOGICAL  UNIT  NUMBER  FOR  OUTPUT  FILE  (20) 

INEXT  -  RECORD  NUMBER  IN  INPUT  FILE  TO  HOLD  NEXT  RECORD  NO. 


r 
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I LAST  -  RECORD  NUMBER  FOR  THE  BEGINNING  OF  THE  LAST  WAVEFORM 
ISIZ  -  MAXIMUM  SIZE  FOR  THE  INPUT  FILE  (  REQUIRED  BY  THE 
RT-11  DEFINE  FILE  STATEMENT  ) 

A  -  SLOPE  FOR  DEVICE  COORDINATE  CONVERSION 
B  -  INTERCEPT  FOR  DEVICE  COORDINATE  CONVERSION 
C  -  CONSTANT  FOR  CALCULATING  CONDUCTIVITIES 
D  -  RATIO  OF  TIME  TO  DEVICE  COORDINATES 


PASS 3  OF  DATA  PROCESSING 


THE  MAIN  PROGRAM  WILL  SET  LOGICAL  UNIT  NUMBER  5  AS  THE  STANDARD 
INPUT  (  TERMINAL  )  AND  6  AS  THE  STANDARD  OUTPUT  (  TERMINAL  ) 

COMMON  DECLEARATION 

COMMON  /DATA/  IDATA(4096)  !  HOLD  INPUT  DATA 

COMMON  /PARM/  IIN, IOUT.INEXT, ILAST, ISIZ,A,B,C,D  !  PARAMETERS 

FOR  ALTITUDE  INTERPOLATION 

COMMON  /ALTITU/  NALT,T( 100) ,Z( 100) 

EQUIVALENCE  PART  OF  THE  HEADER  INFORMATION 

TIMING,  NUMBER  OF  (PAIR)  POINTS,  NUMBER  OF  256  WORD  BLOCKS 

EQUIVALENCE  (ITIME , IDATA( 121 )), (NPTS, IDATA( 125) ) .(NBLK, IDATA< 127) ) 

DATA  ICHAR, JCHAR  /'  ','Y  •/  !  NOTE  UNIX  WILL  PUT  'Y'  IN  THE  LOWER  BYTE 

START 

5  FOR  TERMINAL  INPUT  IS  DEFAULT  FOR  RT-11 

CALL  ASSIGN(6  'TT.-/N')  !  ASSIGN  6  FOR  TERMINAL  OUTPUT 

UNIX  DO  NOT  REQUIRE  THIS,  THIS  IS  DEFAULT 

NWAVE  =  1  !  WAVEFORM  NUMBER 

WRITEC6, 1 )  !  WRITE  OUT  MESSAGE 

FORMAT^ »»*»  PASS3  OF  DATA  PROCESSING  ««') 

INITILIZE  PARAMETERS 
CALL  INIT 

DETERMINE  WHICH  WAVEFORM  TO  START 
WRITE(6,2) 

FORMAT ( 'WHICH  WAVEFORM  TO  START?, 14  FORMAT') 

READ (5. 3)  NSTART 
FORMAT (I 4) 

IF (  NSTART  .LE.  1  )  GOTO  10  !  A  <CR>  IS  INTERPERATE  AS  0 

OTHERWISE  SKIP  WAVEFORMS 

DO  100  I  =  1.NSTART-1  !  SKIP  NSTART  -  1  WAVEFORMS 

READIN  WILL  ftEAD  IN  THE  NEXT  WAVEFORM 
CALL  READIN 

NWAVE  =  NSTART  !  UPDATE  WAVEFORM  NUMBER 

THE  MAIN  LOOP  IS  HERE 

CALL  READIN  !  READ  IN  THE  DESIRED  WAVEFORM 

THIS  IS  WHY  NSTART-1  WAVEFORM  IS  SKIPPED  BECAUSE  THIS  CALL  TO 
READIN  WILL  GET  TO  THE  RIGHT  WAVEFORM 

CALL  LED( ITIME)  !  DISPLAY  THE  TIMING  OF  THIS  WAVEFORM 

WRITE  SOME  INFORMATION  OUT 
WRITE(6, 4)  NWAVE, ITIME, NPTS, NBLK 


FORMAT ( '  WAVEFORM  ' ,  14. 5X ,  ’TIME  =  M5 , 5X , 

'DATA  POINTS  z ’,15. 5X  'BLOCKS  =*.I4) 

IF  THIS  IS  AN  EMPTY  WAVEFORM,  SKIP  TO  THE  NEXT  ONE 

IF(  NPTS  .LE.  10)  GOTO  11  !  UPDATE  COUNTER  AND  READ  AGAIN 

DISPLAY  WAVEFORM  WITH/WITHOUT  ANALOG  CHANNEL 

CALL  DISPRC1, NPTS, 0,200,1)  !  ALL  POINTS.  0-200  HZ,  FLAG  z  1 

CALL  SLOPE  !  GET  BREAK  POINTS  AND  FIT  SLOPE 

SEE  IF  THERE  ARE  ANY  MORE  DATA  FROM  THE  INPUT  FILE 
IF  C  INEXT  .GE.  ISIZ  )  GOTO  1000  !  END 
OTHERWISE  ! 


NWAVE  z  NWAVE  +  1 
WRITEC6.5) 

FORMAT( 'CONTINUE  ?') 

READ (6, 6)  ICHAR 
FORMAT(AI) 

IF (  ICHAR  .EQ.  JCHAR  )  GOTO  10 
OTHERWISE 
WRITE(6,7) 

FORMAT ( *•"••  END  OF  PASS 3  »••••) 
END 


!  UPDATE  WAVEFORM  NUMBER 


!  GET  ANSWER  FROM  THE  TERMINAL 


!  YES,  CONTINUE 


!  END  MESSAGE 


THIS  PROGRAM  MODULE  CONTAINS  MOST  OF  THE  SUBROUTINES 
PASS 3  CALLS 


SUBROUTINE  INIT  INITILIZES  PARAMETERS  AND  OPEN  INPUT  AND 
OUTPUT  FILES 

SUBROUTINE  INIT 
COMMON  DEFINITIONS 
COMMON  /DATA/  IDATA(4096) 

COMMON  /PARM/  IIN. IOUT. INEXT, ILAST, ISIZ, A,B,C,D 
COMMON  /ALT II U/  NALT.TC 100) ,Z( 100) 

INITILIZE  I/O  LOGICAL  UNIT  NUMBERS 
IIN  =  19 
IOUT  z  20 

INPUT  ALTITUDE  INFO 
WRITE(6,5) 

FORMAT( 'INPUT  TIME-ALTITUDE  FILE  NAME',/) 

THIS  IS  A  RT-11  SYSTEM  SUBROUTINE 
CALL  ASSIGN (2 1 , '???.???' ,-1 , 'RDO' ) 

DO  100  I  z  1,100 

READ  IN  ASCII  RADAR  DATA  FILE 

READ(21,10,ENDz101)  T(I),Z(I) 

CONTINUE 
FORMAT (2F8. 2) 

OPEN  I/O  FILES 

GET  THE  NUMBER  OF  RADAR  DATA  POINTS 
NALT  z  1-1 

UNIX  WILL  CALL  "SETFIL"  WHICH  DO  THE  SAME  THING 


OPEN  "PASS2.TMP" 


CALL  ASSIGN(IIN, 'PASS2.TMP' ,9. 'OLD')  I  INPUT  FILE 
CREATE  "PASS3.TMP" ,  ASCII  OUTPUT 

CALL  ASSIGN ( IOUT, 'PASS3.TMP' ,9, 'NEW' )  !  CREATE  OUTPUT  FILE 

A  AND  B  ARE  CONSTANTS  TO  CONVERT  FREQUENCY  TO  DEVICE  COORDINATE 
ON  THE  TEKTRONIX  603 
A  z  2.0475000E1 

B  z  0.0  !  START  WITH  0  -  200  HZ 
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I,  i  i 


WRITE(6, 1 ) 

FORMAT (J  NUMBER  OF  BLOCKS  IN  PASS2.TMP,  14  FORMAT') 

READ(5,2)  N 
FORMAT (14) 

RT-11  DEFINE  FILE  STATEMENT  REQUIRES  THIS 
ISIZ  s  N  I  SAVE  FILE  SIZE 

THIS  THE  IS  RT-11  FORMAT 
SOME  OTHER  SYSTEM  MIGHT  NEED  MODIFICATIONS 

DEFINE  FILE  IIN(  N,  256,  U,  INEXT  )  !  INEXT  WILL  POSITION  TO  THE 

I  NEXT  AVALIABLE  RECORD 

WRITE(6, 3) 

FORMAT ( 1  CALCULATION  CONSTANT? ' ) 

SEE  APPENDIX  B  FOR  FURTHER  DETAIL 

READ<6,4)  C  !  SIGMA  =  C  •  DF/DF 

FORMAT (E 12. 5) 

INEXT  =  1  t  MAKE  SURE  IT  STARTS  READING  FROM 

BLOCK  1 

RETURN 

END 

THIS  ROUTINE  READS  IN  THE  NEXT  WAVEFORM 

INEXT  ALWAYS  POINTS  TO  THE  NEXT  RECORD  AVAILABLE 

SET  I LAST  <-  INEXT  TO  REMEMBER  WHERE  THE  LAST  WAVEFORM  WAS 

SO  THAT  AFTER  THE  CALCULATIONS  THE  INFORMATION  CAN  BE 

WRITTEN  BACK  TO  THE  HEADER 

READ  IN  THE  1 'ST  BLOCK  FIRST  BECAUSE  EACH  WAVEFORM  IS  AT 
LEAST  1  BLOCK  LONG 

FIND  OUT  HOW  LONG  THIS  WAVEFORM  IS  AND  READ  IN  THE  REST 

SUBROUTINE  READIN 
COMMON  DEFINITIONS 
COMMON  /DATA/  IDATA(4096) 

COMMON  /PARM/  IIN.IOUT.INEXT.ILAST.ISIZ.A.B.C.D 
EQUIVALENCE  THE  1 'ST  256  WORDS 
INTEGER  IBUF ( 256 ) . JBUF ( 256 ) 

EQUIVALENCE  ( IBUF ( 1 ) , IDATA ( 1 ) ) , ( NBLK , IDATA (127)) 

SET  ILAST  =  INEXT 

ILAST  =  INEXT 

READ  IN  THE  1'ST  BLOCK 

READUIN' INEXT)  IBUF  I  READ  IN  256  WORD  RECORD 

IF(  NBLK  .EQ.  1  )  RETURN  !  ONLY  1  BLOCK  LONG 


OTHERWISE,  READ  IN  THE  REST  OF  THE  WAVEFORM 


I  lli 

ISTART  =  256 


!  THERE  ARE  256  POINTS  ALREADY 
I  START  AT  THE  257  TH  POINT 
DO  100  I  =  1 , NBLK  -  1  I  DO  THE  REST 

READdIN' INEXT)  JBUF  !  READ  IN  256  WORDS  A  TIME 

COPY  JBUF  INTO  IDATA  AT  THE  RIGHT  BOUNDARY 
DO  101  J  =  1,256 
IDATA (J  +  ISTART)  =  JBUF(J) 

ISTART  =  ISTART  +  256  !  SET  FOR  THE  NEXT  256  WORD  BOUNDARY 

CONTINUE 

THE  WHOLE  WAVEFORM  IS  IN  IDATA(»),  CAN  RETURN  NOW 

RETURN 

END 

THIS  ROUTINE  EXPANDS  THE  WAVEFORM  TO  A  NEW  SET  OF  LIMITS 
AND  DISPLAYS  WITHOUT  THE  ANALOG  DATA 

SUBROUTINE  EXPAND 

COMMON  /PARM/  IIN.IOUT, INEXT, ILAST, ISIZ, A.B.C.D 

DATA  ICHAR, JCHAR  /'  ','Y  '/  J  UNIX  PUT  *Y*  IN  THE  LOWER  BYTE 
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10 

1 

2 

C 

3 


C 

C 


10 


!  A  <CR>  MEANS  NO 


I  GET  THE  LIMITS 


WRITE(6,1) 

FORMAT ('EXPAND  ?') 

READ(5, 2)  ICHAR 
FORMAT(AI) 

IF( ICHAR  .NE.  JCHAR )  RETURN  !  DO  NOT  EXPAND 
OTHERWISE  I 
WRITE(6, 3) 

FORMAT ('LOW-X,  HIGH-X,  LOW-Y,  HIGH-Y  ?•) 

READ(5,4)  ILX, IUX, ILY, IUY 
FORMAT( 41 10) 

CALCULATE  NEW  A  AND  B 

FOR  DEVICE  COORDINATES  EXPANSION 

A  =  4095. /( FLOAT (IUY-ILY)) 

B  =  -A  «  FLOAT (ILY) 

DISPLAY  THE  EXPANDED  WAVEFORM  WITHOUT  THE  ANALOG  CHANNEL 
CALL  DISPR ( ILX , IUX , ILY , IUY , 0 ) 

GOTO  10  !  AGAIN? 

RETURN 

END 

THIS  ROUTINE  CALCULATE  THE  CONDUCTIVITIES  AND 
OUTPUTS  THEM  TO 
THE  OUTPUT  FILE 

ALSO  PUT  THESE  VALUES  BACK  INTO  THE  HEADER  AND  RESETS  ALL 

NECESSARY  INFORMATION 

SUBROUTINE  ENDW 

COMMON  /DATA/  IDATA(4096) 

COMMON  /PARM/  UN. IOUT. INEXT, ILAST, ISIZ , A.B.C.D 
THIS  IS  EQUIVALENCED  TO  THE  HEADER 
INTEGER  IBUF (256) 

EQUIVALENCE  ( IBUF ( 1 ) , IDATA ( 1 ) ) 

EQUIVALENCE ( ITIME , IDATA ( 1 2 1 ) )± ( ALTI , IDATA (41)) 


TIME  =  FLOAT (ITIME) 

PASS  IT  TO  THE  ALTITUDE  INTERPOLATING  ROUTINE 
ALTI  =  ALT (TIME) 

I  =  INEXT 

GO  BACK  TO  THE  LAST  WAVEFORM  AND  WRITE  OUT  THE  HEADER 
WRITE (I IN 'ILAST)  IBUF  !  WRITE  THE  HEADER  BACK 

GO  FORWARD  TO  WHERE  THE  PRESENT  WAVEFORM  IS 
INEXT  =  I 

FIND(IIN'I)  I  POSITION  BACK  TO  THE  NEXT  RECORD 

A  =  2.0475000E1  !  RESET  A  AND  B  TO  0  -  200 

B  =  0.0 

WRITE  OUT  THE  INFORMATION  TO  THE  SYSTEM'S  CONSOLE 
AND  THE  FILE  "PASS3.TMP" 

WRITE ( IOUT , 10 )  ITIHE . ALTI , SIGP  I , SIGN  1 , SIGP2 
, SIGN 2, IDATA ( 5 ) , IDATA (6 ) 

WRITE (6^10)  ITIME , ALTI , SIGP 1 , SIGN  1 , SIGP2 , SIGN2 , IDATA ( 5 ) , 

FORMAT (I 4, 1X,F6.2,4(1X,E12.5),11X,A2,I3) 

RETURN 

END 


r 


THE  FOLLOWING  TWO  ROUTINES  PROVIDE  LED  DISPLAY  AND 
TEKTRONIX  603  DISPLAY 

LED  AND  LEDD  ARE  TWO  ENTRY  POINTS  FOR  DISPLAYING  NUMERIC 
VALUES  ON  THE  LED  DISPLAY 

LED  IS  CALL  FROM  FORTRAN  AND  LEDD  EXPECTS  THE  VALUES  IN  R4  ON  CALLING 


.NLIST 

.MACRO 

MOV 

MOV 

MOV 

MOV 

MOV 

MOV 

.ENDM 

.MACRO 


PUSH 

%0, -(.%(> ) 

*1  -  *6} 

*3, -(*6) 


POP 


MOV 

*6)+ 

MOV 

*6  + 

MOV 

(X6)  + 

MOV 

(%6)  + 

MOV 

X6U 

MOV 

(X6 )  + 

.ENDM 

.’macro 

RETURN 

RTS  *7 

.ENDM 

! MACRO 

CALL  A 

JSR 

X7.A 

.ENDM 

,%3 


LED: 


LEDD: 

START: 


.LIST 
;  START 

. MCALL  . . V2 . . , . REGDEF 
. .V2.. 

.REGDEF 

;  DECLARE  ENTRY  POINTS 
IGLOBL  LED  LEDD 
.TITLE  LED  DISPLAY 
:  LPS  ADDRESSES 
ADBF=170402 
.CSECT  LED IS 


PUSH 
MOV 
BR 

:  THIS 
PUSH 

CLR 


@2(R5),R4 

START 

IS  THE  WAY  FORTRAN  PASSES  PARAMETERS 
;  IN  R4  ALREADY 

LED1 


CLEAR  ALL  6  LED  SEGMENTS 
MOV  #17,e#ADBF 

MOV  #417.§#ADBF 

MOV  #1017,f#ADBF 

MOV  #1417,|#ADBF 

MOV  #2017.«#ADBF 

MOV  #2417,|#ADBF 

;  COUNTER  FOR  6  SEGMENTS 
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1$: 

2$: 


LED1 : 
LED2: 
LED3: 


MOV  #-1,LED2 

INC  LED2 

:  GET  THE  MODULO  TO  CONVERT  FROM  BINARY  TO  DECIMAL 
SUB  #10., R4 

BGE  2$ 

ADD  I10..R4 

MOVB  R4, l£d3 

MOVB  LED1.LED3+1 
LIGHT  UP  1  SEGMENT 


MOV 
INC 
MOV 
TST 
BNE 
POP 
RETURN 
WORD  0 
WORD  0 
WORD  0 
END  LED 


CSECT  ERASVC 


LED3,##ADBF 
LED1 
LED2.R4 
R4 
1$ 


ERASE: 

1$: 


ERASE  AND  SET  UP  THE  STORAGE  SCOPE 

SBTTL  ERASE 
LPSVC  STATUS 
VCST  =  170416 
:  SET  ERASE  BIT 
ERAS  =  10000 

;  Y  MODE,  FAST  INTENSIFY 
VCRDY  =2010 
.GLOBL  ERASE 

;  NAME  CALLED  BY  FORTRAN 

;  ENTRY  POINT 

PUSH 

MOV  #ERAS,«#VCST 
MOV  #VCRDY,€#VCST 

TSTB  €#VCST 

BPL  1$ 

POP 

RETURN 


ROUTINE  TO  INTENSIFY  AN  ADDRESSED  POINT  ON  THE  603 
X-Y  COORDINATES  ARE  PASSED  FROM  FORTRAN 
CSECT  INTE 

SBTTL  INTENSIFY  A  POINT 
POINT 
INTEN 

X  COORDINATE 
Y  COORDINATE 


INTEN: 


2$: 


ENTRY 
GLOBL 
LPSVC 
,  LPSVC  * 

VCX  =  VCST+2 
VCY  =  VCST+4 
PUSH 

;  GET  X  AND  Y 
MOV  |2(R5 

MOV  f4(R5 

:  TEST  FOR  RE 
TSTB  eiVCST 
BPL  2$ 

POP 

RETURN 


VALUE  INTO  THE  LPSVC  REGISTERS 
~#V  CX 
#VCY 
BIT 


.END 

THIS  ROUTINE  DISPLAYS  THE  WAVEFORH  CONTAINED  IN  IDATA(») 
ACCORDING  TO  THE  LIMITS  EXPRESSED  BY  1X1 . 1X2. IY1 , IY2 
IFLAG  IS  USED  TO  INDICATE  WHEATHER  THE  A&ALO&  CHANNEL 
CHANNEL  1  OF  THE  LPS  ADC  IS  TO  BE  DISPLAYED  OR  NOT 

SUBROUTINE  DISPR ( IX 1 . 1X2 . IY 1 , IY2 . IFLAG ) 

COMMON/DATA/  IDATA(4&96) 

COMMON /PARM/  IIN. IOUT . INEXT , ILAST, IS IZ , A , B , C , D 
THIS  IS  TO  HOLD  THE  VERTICAL  INTERVALS 
DIMENSION  IFF(9) 

CURRENTLY  IS  MARKED  AT  EVERY  20  HZ  INTERVAL 
DATA  IFF/  20.40.60.80,100,120,140,160,180/ 

THE  NUMBER  OF  POINTS 

N  =  1X2  -  1X1  +  1 

SKIP  THE  HEADER 

1ST  =  129  +  2  *  (1X1  -  1) 

ERASE  THE  SCOPE 

CALL  ERASE 

DO  100  I  =  1,9 

DRAW  THE  VERTICAL  GRIDS 

IY  =  IFF(I) 

IF (  IY  .LE.  IY1  )  GOTO  100 
IF(  IY  .GE.  IY2  5  GOTO  101 
IY  =  INT(  A  *  FLOAT(IY)  +  B  ) 

DO  102  J  =  1,4096,25 

IX  =  J  -  1 

CALL  INTENUX.IY) 

CONTINUE 

SUM  UP  THE  TOTAL  TIME  FOR  HORIZONTAL  EXPANSION 
TTOTAL  =0.0 
I  =  1ST 

DO  106  J  a  1.N 

TTOTAL  =  TTOTAL  +  COUNT (IDATA(I)) 

1  =  1  +  2 

TTOTAL  =  1.0E-5  *  TTOTAL 

GET  THE  SCALING  FACTOR  IN  THE  X  AXIS 

D  =  4096.0/TTOTAL 

T  =  0.0 

DRAW  OUT  HORIZONTAL  GRIDS 
AND  DATA  POINTS 

MARKING  IS  DONE  ON  EVERY  50  POINTS 

DO  104  I  =  1.50 

CNT  =  COUNT (IDATA( 1ST)) 

F  =  1.0E5/CNT 
T  =  T  +  1 . 0E-5*CNT 
Y  =  A  •  F  +  B 

IF(  Y  .LE.  0.0  .OR.  Y  .GE.  4095.0  )  GOTO  105 
IX  =  INT ( (  D  •  T  ) ) 

IY  =  INT(Y) 

CALL  INTENUX.IY) 

SKIP  THE  ANALOG  POINT 
1ST  =  1ST  +  2 
N  “  N  •  1 

IF(  N  7EQ.  0  )  GOTO  110 
CONTINUE 

THIS  PART  REALLY  DRAWS  OUT  THE  HORIZONTAL  LINE 
DO  108  I  s  1,4096,25 
IY  =  I  -  1 
CALL  INTENUX.IY) 

GOTO  107 

EXPAND  LATER  TO  DRAW  OUT  THE  ANALOG  CHANNEL  ALSO 
RETURN  !  EXPAND  LATER 
END 
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THIS  FUNCTION  CONVERT  THE  COUNT  INTO  A  REAL  NUMBER 
SINCE  THE  COUNT  USES  FULL  16  BITS,  IF  THE  COUNT 
EXCEED  15  BITS,  IT  WILL  BE  VIEWED  AS  A  NEGATIVE  INTEGER 
BY  FORTRAN,  THEREFORE  THE  PROGRAM  HAS  TO  CHECK  FOR  THAT 
NEGATE  THE  NUMBER  AND  ADD  32767  TO  IT 
REAL  FUNCTION  COUNT (ICNT) 

IF(  ICNT  .LE.  0  )  GOTO  1 
COUNT  =  FLOAT (ICNT) 

RETURN 

COUNT  =  FLOAT (-ICNT)  +  32767.0 

RETURN 

END 


THIS  ROUTINE  EXPONENTIALLY  INTERPOLATES  THE  ALTITUDE 
FROM  THE  TIME  GIVEN 

ALTITUDE  DATA  IS  KEPT  IN  COMMON  FROM  WHICH  IS  READ  IN 
FROM  THE  RADAR  DATA  FILE 

THE  NUMBER  OF  RADAR  DATA  POINTS  IS  THERE  ALSO 
REAL  FUNCTION  ALT (TIME) 

COMMON  /ALTITU/  NALT,T( 100) ,Z( 100) 

IF (  TIME  .LT.  T(1))  GOTO  10 

DO  100  I  =  2 , NALT 

FIND  THE  RIGAT  VALUES 

IF(  T ( I )  .GT.  TIME  )  GOTO  11 

CONTINUE 

T1  =  T ( 1—1 ) 

T2  =  T(I) 

Z1  =  Z(I-1) 

Z2  =  Z(I) 

LINEAR  INTERPOLATE  NOW 

ALT  =  Z1  +  (Z2-Z1)/(T2-T1)*(TIME-T1) 

RETURN 
WRITE(6. 1 ) 

FORMAT ('ERROR  IN  ALTITUDE') 

ALT  r  -1.0 

RETURN 

END 


THIS  ROUTINE  PERFORMS  WEIGHTED  LEAST  SQUARE  STRAIGHT  LINE  FIT 

FOR  THE  SLOPE  OF  THE  WAVEFORM 

BREAK  POINTS  ARE  INPUTTED  FROM  THE  TERMINAL 

FIRST  TRAIL  OF  THE  FIT  ASSIGNS  EQUAL  WEIGHTS(I.O)  TO  ALL  THE  POINTS 
WITHIN  THE  INTERVAL 

SECOND  TRIAL  ASSIGNS  WEIGHTS  BETWEEN  0.0  AND  1.0  TO  THE  POINTS 
ACCORDING  TO  THEIR  RESIDUE  BETWEEN  THE  FIRST  FIT  AND  THE  ACTUAL 
VALUE 

THIRD  TRIAL  WILL  DO  THE  SAME  THING  TO  FINER  TUNE  THE  LINE 
THE  FITTED  LINE  IS  THEN  DISPLAYED  ON  THE  SCOPE 

CONDUCTIVITY  VALUES  ARE  THEN  CALCULATED  AND  INSERTED  INTO  THE  HEADER 
THE  PROCESS  IS  REPEATED  FOR  NEGATIVE  CONDUCTIVITY  AND  MORE 
IF  THERE  WERE  MORE  THAN  ONE  SET  OF  CONDUCTIVITY  VALUES 


DATA  BASE  DEFINATION 


SUBROUTINE  SLOPE 
COMMON/DATA/  IDATA(N096) 

COMMON  /PARM/  IIN, IOUf ,INEXT,ILAST,ISIZ,A,B,C,D 
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REAL  SIG<8,4)  !  THIS  IS  1 

LOGICAL  FUG  !  FUG  TO  : 

EQUIVALENCE  IT  TO  THE  HEADER 
EQUIVALENCE  (SIG( 1 , 1 ) ,IDATA(49)) 

DATA  ICHAR, JCHAR/  *  ',»Y  V  !  UNIX  PUT 

FLAG  =  .FALSE.  !  FALSE  =  ] 

SEE  IF  EXPANSION  IS  NEEDED 

CALL  EXPAND 

GET  THE  BREAK  POINTS 

WRITE(6. 3) 

FORMAT ('LOW-X,  HIGH-X  ?') 

READC5.4)  1X1 ,1X2 
F0RMAT(2I 10) 

IF (  1X1  .EQ.  0  )  GOTO  111  !  A  <CR>  = 

FOR  NEGATIVE  COND.  IF  YES.  SIGMA+  r  SIGMA - 
IS  THE  LIMITS  PUT  IN  CORRECT 
IF(  1X2  .GT.  1X1  )  GOTO  1 
WRITE (6, 2) 

FORMAT (  •  ERROR  IN  LIMIT  *) 

GOTO  150 


I  THIS  IS  THE  SIGMA  BLOCK  FROM  HEADER 
!  FUG  TO  INDICATE  POS  OR  NEG 


f  UNIX  PUT  IN  THE  LOWER  BYTE 
I  FALSE  =  POSITIVE;  TRUE  s  NEGATIVE 


t  A  <CR>  =  0,  SEE  IF  NOW  IS  LOOKING 


NUMBER  OF  POINTS 
N  s  1X2  -  1X1  +  1 
JSTART  =  129  +  2»(IX1 
INDEX  =  1 
II  =  1 
EPS  =  1.0 
TS  =  0.0 


I  TAKE  POINTS  FROM  1X1  -  1X2  INCLUSIVE 
1 )  I  SKIP  THE  HEADER  AND  GO  TO  THE  RIGHT  START 
!  WHICH  PASS 

!  FOR  FILLING  IN  SIGMA  VALUES 
!  TOLERANCE  FOR  RESIDUE 
f  STARTING  TIME 


HAVE  TO  SKIP  THE  DATA  BETWEEN  POINT  1  AND  1X1  AND  CALCULATE  THE  TIME 
ELAPSED 

DO  110  I  s  129. JSTART, 2  t  POINTS  ARE  IN  PAIRS 
SUM  UP  THE  COUNTS 


DO  110  I  =  129. JSTART, 2  t 
SUM  UP  THE  COUNTS 
TS  =  TS  +  COUNT (IDATA(I)) 
TS  =  1.0E-5  *  TS  ! 

T  a  TS  ! 

STARTING  INDEX 
J  =  JSTART 

FREQI  =  0.0  ! 

WI  =  0.0  I 

TF  =  0.0  ! 

TI  =  0.0  l 

TI2  =  0.0  I 

DO  100  I  =  1.N  IN  POINT 
SUM  UP  THE  COUNTS 
CNT  r  COUNT(IDATA(J)) 

F  =  1.0E5/CNT  I 

T  =  T  +  1.0E-5  •  CNT  ! 


!  CONVERT  TO  TIME  ELIPSED 
I  T  STARTS  FROM  TS 


!  SUM  FI 
!  SUM  WI 
I  SUM  TI  *  FI 
I  SUM  TI 
I  SUM  TI  ••  2 
I  N  POINTS  TOTAL  ,  ST 


START  FROM  THE  J  TH 


I  CONV  TO  FREQUENCY 
I  TIME 


GOTO  (101,102,102,103),  INDEX  I  WHICH  PASS  IT  IS 
W  s  1.0  I  FIRST  TIME  WEIGHT  r  1  FOR  ALL  POINTS 

GOTO  105 

R  =  ABS(P*T  ♦  Q  -  F)  f  RESIDUE 
IF (  R  .LE.  EPS  )  GOTO  107 
CALCULATE  NEW  WEIGHT 
W  “  EPS/R 

GOTO  105  I  CAL  WEIGHT 

W  =  1.0 

SUM  UP  ALL  INFORMATION 

NOTE  THAT  THEY  HAVE  ALL  BEEN  MODIFIED  BY  WEIGHTS 
THE  1 'ST  PASS  ALL  WEIGHTS  ARE  ASSIGNED  1.0  ^ 

FOR  SUCCESSIVE  PASSES  THE  WEIGHTS  ARE  MOfilFIED 

ACCORDING  TO  THEIR  RESIDUES 

WI  s  WI  +  W  I  SUM  WEIGHT 

TI  *  TI  +  W  *  T  I  SUM  TI 

TI2  s  TI2  ♦  W  •  T  ••  2  I  SUM  TI  ••  2 
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+  W 


•  F 

•  F 


100 

C 


c 

103 

C 


120 

121 

C 

C 

200 


130 

201 

9 


112 

C 

6 

7 

C 


140 

141 


SUM  T  *  F 

SUM  FREQUENCY 

DATA  ARE  STORED  IN  PAIRS 

I  DETERMINT 


I  NEXT  PASS 


TF  =  TF  +  W  » 

FREQI  =  FREQI 
J  =  J  +  2 
CONTINUE 

DET  =  WI  *  TI2  -  TI  ••  2 
FIND  SLOPE  AND  INTERCEPT 
P  =  (WI  »  TF  -  TI  *  FREQI )/DET 
Q  s  (TI2  •  FREQI  -  TI  •  TF ) /DET 
INDEX  =  INDEX  +  1 
GOTO  5 

DISPLAY  THE  FINAL  FITTED  LINE 
IF (  FLAG)  GOTO  200 

RESTRICT  TO  200  POINTS  FOR  LINE  GENERATION 
DT  =  20.48  /  D 
T  =  0.0 

DO  120  I  =  1,200 
F  =  P*T  -*■  Q 
Y  =  A*F  +  B 

IF(  Y  .LE.  0.0  )  GOTO  120 
IF (  Y  .GE.  4095.0)  GOTO  121 
IY  =  INT(Y) 

IX  =  INT(D*T) 

CALL  INTEN(IX.IY) 

T  =  T  +  DT 
GOTO  201 

CHOOSE  BETWEEN  WHICH  DIRECTION  TO  GENERATE  THE  LINE  FOR 
BEST  EFFECT 
TTOTAL  =  4096.0/D 
DO  130  I  =  1,201 
F  =  FLOAT (1-1) 

T  =  (F-Q)/P 

IF (  T  .LE.  0.0)  GOTO  130 
“  T  .GE.  TTOTAL)  GOTO  201 
A  »  F  +  B 
INT  (  T  »  D  ) 

IY  =  INT(Y) 

CALL  INTEN(IX.IY) 

CONTINUE 
WRITE(6. 9) 

FORMAT ('TRY  AGAIN?') 

READ(5,7)  ICHAR 

IF (  ICHAR  .EQ.  JCHAR  )  GOTO  150 

CALCULATE  THE  RIGHT  SIGMA  VALUES  AND  PUT  THEM  IN  THE  RIGHT  PLACE 


IF( 

Y  = 
IX  = 


! 


IF(  FLAG  )  GOTO  112 
SIG(I.II)  =  P  »  C  ! 

FLAG  =  .NOT.  FLAG  ! 

GOTO  150  ! 

SIG(5, II )  =  P  *  C  I 

FLAG  =  .NOT.  FLAG  ! 

II  =  II  +  1  I 

AT  MOST  4  SET  OF  SIGMA  VALUES  ARE  ALLOWED 
IF(  II  .GT.  4)  GOTO  141 
WRITE(6,6) 

FORMAT ('ANOTHER  SET?') 

READ(5,7)  ICHAR 
FORMAT(AI) 

IF(  ICHAR  .EQ.  JCHAR  )  GOTO  150 
OTHERWISE,  FILL  IN  THE  REST  WITH  OS 
DO  140  I  =  II. 4 
SIG( 1,1)  =  0.6 
SIG(5,I)  =  0.0 
CONTINUE 
CALL  ENDW 
RETURN 


NEGATIVE 
SIGMA  PLUS 
UPDATE  FUG 
TRY  FOR  SIGMA  MINUS 
SIGMA  MINUS 
UPDATE  FUG 

IN  CASE  THERE  IS  MORE  THAN  1  SIGMA  VALUE 


non  o  oi\>  -»  o  ooo  o  rmnnnrjononnonoooooooooooonoooooonooooo 
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111  IF(  FUG  )  GOTO  112 
RETURN 
END 

PROGRAM  TO  EXTRACT  DATA  FROM  HEADER 
INTO  BASIC-PLUS  VIRTUAL  ARRAY  (DISK  FILE)  FORMAT 
FOR  PLOTTING  ON  THE  PDP  11/45 

VIRTUAL  ARRAY  IS  ORGANIZED  AS  FOLLOW; 

IN  BASIC  — 

OPEN  "XXXXX"  AS  FILE  #? 

DIM  #?,  X(300),Y(300> 

REM  X(0)  CONTAINS  THE  NUMBER  OF  COORDINATES 
REM  Y(0)  IS  NOT  USED 

REM  SINCE  BASIC  STARTS  ARRAY  AT  OTH  ELEMENT 
REM  SO  THERE  ARE  ACTUALLY  301  PAIRS  OF 
REM  POINTS  IN  THE  FILE 

REM  THIS  NUMBER  MUST  BE  MAINTAINED  FOR  PROPER 
REM  DATA  ALIGNMENT 


IN  FORTRAN  — 

DISK  FILE  CAN  BE  ADDRESSED  BY  BINARY  I/O 
AND  DEFINE  FILE  STATEMENT 

CALL  ASSIGN (??????????????)  TO  ASSIGN 
LOGICAL  UNIT  NUMBER  FOR  FILE(S) 

DEFINE  FILE  ?(602,4,U,II) 

602  RECORDS  (301  XS  +  301  YS) 

EACH  4  WORDS  LONG  (BASIC  PUTS  EVERYTHING  IN 
DOUBLE  PRECISION) 

READ  (?'IREC)  FOR  READ 
WRITE(? ' IREC)  FOR  WRITE 
FIND(? ' IREC)  FOR  POSITIONING 

SEE  PDP  1 1/RT-1 1  FORTRAN  MANUAL  FOR  REF. 


COMMON  DECLARATION 
COMMON  /DATA/  IDATA(4096) 

COMMON/  PARM/  IIN.IOUT.  INEXT, ILAST, ISIZ, A.B.C.D 
EQUIVALENCE  ( NPTS , IDATA ( 1 25 ) ) , ( ALT , iDATAf  4 1 )  J,( SIGPOS , 
IDATA(49) ) , (SIGNEG, IDATA(57) ) 

THEY  ARE  USED  TO  HOLD  ALTITUDE,  SIG+,  SIG- 
REAL*8  Z.X.Y 

LOGICAL  UNIT  NUMBER  6  IS  INPUT  TERMINAL 
RT-11  DEFAULT  TERMINAL  TO  7 
CALL  ASSIGN (6, 'TT:/N') 

WRITE(6, 1) 

FORMAT ( ' EXTRACT  DATA  FROM  HEADER') 

WRITE(6.2) 

FORMAT (’INPUT  DATA  FILE  NAME  ?',/) 

THE  PASS 2  FILE  NAME 

CALL  ASSIGN ( 19, ' ??????.??? ' ,-1 , 'RDO' ) 

IIN  =  19 

WRITE(6,3) 

THE  3  OUTPUT  FILE  NAMES 
SIGMA+,  SIGMA-,  SIGMAt-  r  SIGMA- 


nno  oo  on  no  -»ooo  non  o  oo 


123 


3 
C 

4 
C 

5 
C 

C 

6 

11 


FORMATC OUTPUT  FILE  1  ~  POS  /) 

CALL  ASSIGN (20, '??????.???' ,-1, JNEW') 

WRITE(6.4) 

FORMAT ('OUTPUT  FILE  2  —  NEG  ?', /) 

CALL  ASSIGN (21 , '??????.???* ,-1 , 'NEW' ) 

WRITE(6,5) 

FORMAT ( ’OUTPUT  FILE  3  —  COM  ?', /) 

CALL  ASSIGN(22, '??????.???' ,-1 , *NEW' ) 

WRITE (6, 6) 

THIS  IS  REQUIRED  BY  RT-1 1  DEFINE  FILE  STM 
FORMAT ('INPUT  FILE  SIZE  ?') 

'll)  ISIZ 


DEFINE  FILE  19(N. 256, U, INEXT) 

DEFINE  FILE  20(602, 4, U, II) 

DEFINE  FILE  21 (602,4, U, 12) 

DEFINE  FILE  22(602, 4, U, 13) 

INEXT  =  1 

SINCE  FORTRAN  STARTS  ARRAY  AT  1 
THEREFORE  BASIC  ARRAY ( 1 )  =  FORTRAN  ARRAY(2) 

IPOS  =  2 
INEG  =  2 
ICOM  =  2 

START  READING 

00  CALL  READIN 

IF (INEXT  .GE.  ISIZ)  GOTO  1000 
IF(NPTS  .LE.  10)  GOTO  100 
IF (  ALT  .LE.  0.0)  GOTO  100 
IF(SIGPOS  .EQ.  SIGNEG)  GOTO  101 

CONVERT  TO  DOUBLE  PRECISION 
Z  =  DBLE ( ALT ) 

X  =  DBLE(SIGPOS) 

Y  =  DBLE (SIGNEG) 

WRITE  THEM  OUT  TO  DISK 
WRITE (20 ' IPOS)  X 
WRITE ( 20 ’IPOS+301)  Z 
WRITE (21 'INEG)  Y 
WRITE(21 'INEG+301 )  Z 

INCREMENT  THE  COUNTER 
IPOS  =  IPOS  +  1 
INEG  =  INEG  +  1 
GOTO  100 


IF  THEY  WERE  EQUAL 
101  Z  =  DBLE (ALT) 

X  =  DBLE(SIGPOS) 

C 

WRITE (22 'ICOM)  X 
WR ITE ( 22  * ICOM+30 1 )  Z 
ICOM  =  ICOM  +  1 
C 


OOOf) 


GOTO  100 


i 

1 


1000 


WRITE  OUT  THE  NUMBER  OF  POINTS 


WRITE 

WRITE 


WRITE 


END 


DBLE (FLOAT ( 
DBLE ( FLOAT c 
DBLE (FLOAT ( 


IPOS-2 

INEG-2 


I COM-2 


ENt 

DATE 
FILMEI 

o-e 

DTIC 


